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ON THE PROBLEM OF THREE BODIES IN A PLANE 


by 
Victor G. Szebehely* 


INTRODUCTION 


In this paper the 12th order system of differential equations of 
the planar three bodies problem is reduced to one of the 4th order by 
introducing generalized coordinates. 

The system of differential equations of the motion for the planar 
problem is one of the 12th order, because the positions of the three 
bolies require 6 coordinates and the kinematic equations are of the 
second order. This system is reducible from the 12th to the &th order 
by using four rntegrals regarding the motion of the center of gravity. 
By means of a contact transformation, this can be reduced to the 6th 
order and using the integral of energy ad eliminating the time, a system 
of 4th order can be obtained as the final result. 

In the present paper, a new set of generalized coordinates is in- 
troduced by means of which the final form of the 4th order system 
can be obtained without a contact transformation. The special case for 
which the masses of the three bodies are equal and the attracting 
forces are proportional to the cubes of the distances is investigated 
in detail and two sets of particular solutions are obtained. One of 
these 1s represented by trochoidal curves, and the other set splits 
into two, each of which can be considered as a generalization of the 
Lagrangian collinear solution. 


GENERALIZED COORDINATES 


In two dimensions the positions of the three bodies are described 
by 6 generalized coordinates. The following coordinate systems are 
introduced: 

(1) the (€,7) system which is fixéd in the plane of motion, 

(2) the (x,y) system which has its ormgin at the center of gravity 
of the three bodies, the coordinates of which are * and Np The x and 
y axes coincide with the principal axes of inertia of the three bodies. 
The (x,y) system is moving wth a translational velocity Po (E+ No) 


and it rotates around its origin with the angular velocity y¥, where 


¥ is the angle from the € to the x axis. 


The motion of the center of gravity is described by £, = 7, = 0. 


Now that the motion of the centroid of the system has been determined, 


* Formerly Associate Professor, Virginia Polytechnic Institute. The views 
expressed in this paper are those of the author and are not to be con- 
strued as those of the Navy Department. 


59 


C= + 


60 MATHEMATICS MAGAZINE (Nov. -Dec, 


the following four generalized coordinates are introduced: ' 
(1) the first radius of gyration, 


t=! 


ls 2 
q, = 7 mx, where M = 2 m 


(2) the second radius of gyration, 


1 :.., , 
q> = [Lz ay, » Slgn (q,q,) = Sign Misi 


(3) the generalized angular coordinate, 


yi [= mix 
= arc tan — = Oy, 
x; . 


2 
2 my; 
where the constants 0; satisfy the equation, tan (0, - 0.) 


Mm, 
(* k and (4) the angle of rotation, w, which has already been in- 


mm 


t roduce_]. 
The (x,, 


coordinates as follows: 


i od lel 
cos(* + 0: ); v5 ° = - sin(x +0,)... (1) 


t 


y,;) coordinates are related to the q,, 94>, * generalized 


(43,8 * 1,2,3) 


EQUATICNS CF MOTION 
The use of the generalized coordinates leads 
pression of the kinetic energy of the system: 


to the following ex- 


aim ( 2) 


y . . . . . . . 
- =t€ + ne + (q° +. q,)(%? + w*) + 49,9,*¥ + q, + <a 

The potential energy has the form: 
( 3) 


V= 


n 
2 
| (q?sin?( +o, )+ q5cos’(* + Oo, BM s<'s 


Here K is a constant and the forces are proportional to the (n - 1)st. 
power of the distances. Formlas (2) and (3) are obtained by el ementary, 


straight forward calcul ations. 
The Lagrangian equations of motion are 


'Radau, R. Comptes Rendus, 68, 1465.1468. 
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d (OL aL 
raed" ° 
dt q; 0q; 


where L = T-Vis the kinetic potential, 
q,, the generalized position coordinates (in our case q,, oe 


* and w) 
qi the generalized velocity coordinates (in our case Sb Ses 
x andy). 
Due to the fact that Vis independent of the generalized velocity 
coordinates, the equations of motion become 


d oT | oT eV 
dt (3q,) 


= ( 4) 


ae 


q,, this is q; - (%* + Y*)q, - 2q,x = 


2 


- (x* + v7 )q, - Qq xv = 





for q; [(q, + q,)% + 29,957] 


for 9; [(qi + a5) + 2q,4,%] = (4) 


t 


Equation (6) is the result of the fact that Y is an ignorable (cyclic) 
coorhinate. From equation (4) can be obtained the integral 


(qi + qa)b + 2q,q,% = 0 


wach corresponds to the integral of angular momentum. Using in addition 
to this, the integral of energy, the equations of motion can be reduced 
to the 4th. order. 


SPECIAL CASES 


Using the system of equations \5), several wel] known particular 
solutions can be obtained such as Lagrange’s collinear solution, 
Lagrange’s equi di st ant solution, 0. Pylarinos’ ' homothetic motion, 
D. Sokolov’ s* particular solution, etc. 

It has already been mentioned that the coordinate y is cyclic. 
The * coordinate is cyclic as far as the kinetic energy is concerned. 


'o. Pylarinos: “Ueber die Lagrangeschen Faelle in verallgemeinerten 
Dreikoerperproblem.” Mathematische Zeitschrift, Vol. 47, pp. 351-372. 
*s. Sokolov: ‘Sur un nouveau cas d’intégrabilité dans le probleme rec- 
tiligne de trois corps.” C. R. (Doklady) Acad. Sci. URSS (N.S.) Vol. 46, 
pp. 95-98. 
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It can be shown that it is also cyclic with respect to the potential 
energy if the masses are equal and the forces vary as the cubes of 
the distances. 

For this case 


on - 27 
x; = V2 q,cos(x + se 0: 9%, * V2 qosin(* + 7%: 


T = 3m/2[& + 72 + (qt + gh)Qi? + ¥?) + aqyq.ed + g? + 92) (7) 


‘oO 
27 Km? 
dle (3q,* + 2q,7*q,7 + 39,4) (8) 
Since the x and the W coordinates are cyclic, 
oT oT 
oe: and 3 = const. 


The last formula for the kinetic energy suggests the introduction of 
four new generalized coordinates: 


P, = L 1. * & 
1 1 | (9) 
?, = b, = W-x 


In terms of these new generalized coordinates the kinetic energy, the 
potential energy and the equations of motion are 
*2 he 


"4 Se a t 
. + Pi Ps + PoP, + P, . Po] 


2 
27 Km 
—— (% + 497%) + 92) 


-18 Kmp,(2% + ney 


er 2 
“18 Kmd, (2, + 207) 
24, 
ad =) d 
PaPo = C2 
Eliminating D, and D, from (12) by means of (13) and using convenient 


units (18 Km = 1) the final form of the fourth order system of di ffer- 
ential equations is obtained: 
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d,-— + o(g + 292) = 0 


2 
$ - 24 0.(g¢ + 99%) = 0 
2 3 a 2? 
2 


4 


PARTICULAR SOLUTIONS 


(1) The above equations are satisfied with constant ?, and ?, for 
which case 
+ p 2n -, - 27 
co s(x + — i), Yi -” a sin(x + — L) (15) 
Ir 3 3 
The positions of the particles in the fixed (£,7) system are obtained 
by considering the translation and rotation of the (x,y) coordinate 
system with respect to the (¢,7) system. 
€. = a,(t - t,) + bg + x,cos Y- y. sin Y) 
>1 € fc) é i i 
\ (14) 
; a(t a b. + x,sin Y + Y.cos y) 
It can be assumed without any loss of generality that the initial 
position of the center of gravity coincides with the origin of the 
fixed coordinate system, in which case 


,e kh #6 


and that the velocity of the center of gravity is zero, in which case 


ae s a, = 0), 
Introducing into formulas (14) the expressions obtained in (15) for 
x, and y; gives 


— D 
vid cos(P, + i) + —2 cos(P, ~ 
/2 v2 
1 : 27 Po : 2 , 
n — sin(P? + —1t) + — sin(P ei 


l /2 3 3 J? 4 


Making use of the fact that 9, and ~, are constants and choosing proper 


initial values, the complex position vector for the j-th particle becomes 


y - + in. = Aet (at*+27/3)) ‘ Bet ( &t-27/3)) 
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1. 2.3. 
0 A2 + 2B?) 
= 2(B? + 2A?) 


Equation (18) represents a trochoidal curve. 
Two theorems can be established regarding this solution, 
3 V3 
(a) The area of the triangle formed by the three bodies is > UN 


= const. 


(b) The sum of the squares of the sides of the triangle is 9( 45 + q5) 


= const. 

According to (E) the distance between the bodies is limited, and 
since in addition by (a) the area is constant (and except for special 
cases different from zero), it follows that the bodies cannot collide. 
Therefore, the motion by Laplace’s definition is stable. Theorems (a) 
and (b) also suggest the use of another set of generalized coordinates, 
q,°, and gq + q°. 

A further analysis of equations (18) and (19) results in a one- 
parametric system of epi-and hypocycloidal curves as orbits. Special 
cases correspond to the Lagrangian collinear and equidistant solutions. 
A special case of the epicycloidal solution is obtained when the three 
bodies lie in fixed positions on a straight line, which is rotating 
with constant angular velocity. The degenerate case of the hypocycl oidal 
solution is a fixed straight line, on which the bodies are moving. 
The equidistant solution can be obtained by letting the hypo-and epi- 
cycloidal orbits coincide. In this case the equilateral triangle formed 
by the three bodies rotates uniformly. The above statements will be 
veri fied in another paper. 

(2) Two further more general collinear solutions can be obtained 
by putting q, = 0. Then 2, = %, = @ and equations (20) reduce to 


DP + 39% = c? (20) 


(21) 


The solution of (20) can be obtained in terms of an elliptic integral: 


P tnt F(o) 

t= ft + os > = , 
‘ [ [(3/204 + 2 + e*/¢ yg" - 3/2¢ - ce] 

0 





whose inverse is ?= f(t). 
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Since q, = 0, the three bodies are located on the rotating x axis and 


as on 
x, = V2 DP cos(x + — t), y; = 0 


In consequence of equation (21) 


Analyzing cases (a) and (b) two new generalized collinear solutions 
can be obtained. 


For case (a), according to equations (13) it follows that 


2 . . _ 
al (P, + +4) = 0) 


[If the trivial solution 2= 0 is disregarded, then Y = 0 and therefore 


const. A repeated use of equations (13) mgves 


ct dt 
f(t) 


ie) 


The solution in the moving coordinate system is 
2 f(t)cosl[ i 
x. = Vv f( t)coste ———— + 
t 1 ‘ f2(t) 
0 
and in the fixed system, 


-t 
/2 f(t)cos y% cos(e, | 


J, 


t 


ie) 





/2 f(t) sin Yo cosle,| 


to 


Introducing the polar coordinates 


and a; 


equations (22) become 
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t dt Qn 
p. = 2 f(t) lal at Sal 


and a, = 


In the above solution the three bodies are moving along a fixed 
straight line which makes an angle ¥, with the € axis, ad the motion 
is described by formula (23). 

In case (b) the solution is 


2 
p, = V2 f(t)cos(x, + =i) 





[ dt 

a.=c —— 

t 1 ; f2(t) 
0 


In this solution the bodies are moving on a rotating straight line 
and are governed by the f(t) function. The variable angular velocity 
of the rotating line is given by c,/f?(t). 

Case (b) can be considered a very general case of collinear motion. 

Regarding the generalization of the above expanded method, the three 
dimensional four bodies problem can be mentioned. When the above in- 
troduced Lagrangian coordinates are used in a more general sense, the 
the 24th. order system of differential equations of the three dimensional 
four bodies problem is reducible to the 12th. order. The reduction is 
again symmetric and does not require any contact transformations. The 
three dimensional three bodies problem, being a degenerate case of the 
four bodies problem, can be of course investigated also with the above 
given method, however no complete symmetry can be expected. 


David Taylor Model Basin 
Navy Department 





THEORY OF BUDGETS BASED ON 
PARABOLIC ENGEL CURVES* 


G. A. Baker 


1. Introduction. In econometrics it is usual to assume that the 
quantities of commodities purchased depen’! only upon the total income. 
Thus if x 1s the total income then the expenditure for the ith item is 
represented by x(x). Suppose that there are n items. Now if x is re- 
garded as a parameter then the functions x, will vary as x assumes 
different values and will trace a arve in the n-dimensional quantity 
space (x,, 5, ... x,). This curve is called an Engel Curve after 
Fmest Engel, a German economist. The Engel Curves can be determined 
empirically by observing the consumption of individuals belonging to 
different income classes. 

It has been generally assumed that the expenditure on a single 
item expressed as a function of total income is linear [1], [2], [3]. 
Allen and Powley recognize that the range must be much restricted for 
the linear relationship to hold. Walc simply assumes that the linear 
relationship holds in some sufficiently restricted region but does not 
investigate the length of the interval in which such an assumption is 
valid. 

Davis discusses only the case for linear curves but suggests an 
equation of the form 
(1.1) xz, = a- be 4* 
where x, 1s the expenditure on a particular item; x is total income; 
a, b, and u are positive numbers with 6 < a, as being suitable over a 
more extended range. The difficulties with this suggestion are that 
the curve to represent the expenditure on the combination of two items 
is not the combination of the representations of the two separate items 
and that an asymptotic limit 1s assumed which may wel! be questioned. 
It is true that only a limited amount of food can be consumed but the 
money value of this food may increase indefinitely. 

It must be realized that categories such as food, clothing, etc. 
are more or less arbitrary and in reality include many more or less 
related but distinct items. Thus in fitting a curve to expenditures on 
food, we are really finding a curve that is the resultant of the com- 
bination of many more particularized expenditure curves. 


* Presented at the fourth annual meeting of the Northern California 
Section of the American Mathematical Association, January 31, 1942, 
University of California at Berkeley. 
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It 1s very inportant that the mathematical models of the economic 
world be as realistic and as extended in scope as possible. The purpose 
of this paper is to give a more realistic and extended basis for 
the theory of budgets. 

Lately a great deal of budgetary data has |ecome available as a 
result of the Consumer Purchases Study carried on as a joint project 
by the Bureau of Home Economics and Bureav of Labor Statistics with 
the cooperation of the National Resources Planning Poard, the Works 
Progress Adninistration, and the Central Statistica! Poard. 

An extensive examination of these data indicate that a least-square 
parabola of secon! degree will adequately represent the functional re- 
lation between expenditure on an item and total income or expenditure 
on all items for a very extended range of the income variable. Para- 
bolas have the advantage that the representations for the separate ites 
can be added for the representation of the total of any subset of 
items or the total of all of the separate items. Jf the functional 
relation of expenditure on an item to total expenditure is parabolic 
then the Engel Curve as defined by Wald, loc. cit., page 145, equations 
(2) is paratolic for two expenditures. If only two expendiwres are con- 
sidered and the Engel Curve 1s paralolic and the utility or andicator 
function quadratic, then the utility function is determned except for 
a constant multiplier and an additive constant. The degree of deter- 
minateness is much less in the case of n items since the utility functio 
then involves an arbitrary positive definite quadratic form of n -2 
variables. 

2. Engel curves and utility function. Let us consider two expen- 
ditures, x, = expenditure on savings and x, = expenditure on all other 
items. Then x = x, + x, is the total income or expenditure. For a par- 
ticular set of data taf, given in detail in table 1 we have 


(2.1) x, .0000545 x* + .0433 x - 61.3 


x, = -.0000545 x* + .95467 x + 61.3 


Eliminating x between the equations of (2.1) we obtain 


(2.2) (x, + x)? + 7.9450 x 107(x, + x) - 1.8349 x10*x, 


-1.1248 x 10° = 0 


as the Engel Curve, which is a parabola. It should be noted that the 
Engel Curve can be parabolic only if the prices paid for essentiall) 
the same items differ at different income levels. 

In Davis’ (loc. cit. pp. 166-168) notation the parabolic Engel 
Curve is given by 
}?2 


(2.3) 20,(l.x, + mx.) + rs K(o,(l,x, + mx.) + 
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TABLE 1 


Expenditure Patterns of Southeast Village Families Containing only 
Husband and Wi fe. 


Expenditure Expeiditure on 
Average No. of on savings all other items 
Total Incove families Actual Compu ted Actual Computed 








40 4 18 -3 -34 407 
57 54 -4 -18 585 
788 63 -3 4 7 9] 
1047 55 58 44 Gag 
1303 49 A4 Rg 1 239 
1508 67 117 1 28 1391 
1743 38 195 180 1548 
20469 50 268 262 1201 
2511 19 333 39] 2178 
3087 17 453 592 2434 

é 1212 1142 3143 

T 3092 3114 4152 


if the utalaty function is 


2 2 
(2.4) U(x,, x,) = Ax + 2Bx,x, + Cx 


, , * 2Dx, + 2Exz, + F. 
Take 1, = 1, = 1 and we find that 
= -4.7525 x 107° 
-4.5871 x 10° 
1.9221 x 10’ 
-1.8683 x 10’ 


The values of the coefficients of (2.4) are then determined as 
- 3.4403 x 10° D = 1.3450 x 10° 
1.1468 x 10° E = 9.4760 x 10° 


-3.4403 x 10° F = arbitrary 
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Hence a utility function is 
(2.7) U(x,,%,) = -3.4403x10%x,? + 2.2936 x 10%x, x, 
3.4403 x10°x,* + 2.6900 x10°x, + 1.8952x10’x, + F. 


U,, = ~6.8806 x10* <0 and 


U, ,U,>| |-6-8806 x 10% - 2.2936 x 10° 
- > 0 

3 = 3 . s 
Ler -2. 2936 x 10 6.8806 x 10 : 
Hence U is a convex quadratic function and is suitable as a utility mn 
function. me 
The indifference lines are arcs of ellipses. The axes of these by 
ellipses make a angle of 45° wth the x, and x, axes and the center of = 


these ellipses is the vertex of the parabola (2.2). 

3. The n dimensional case. Professor Charlies B. Mprrey, Jr. of the 
Department of Mathematics of the University of California has considered 
the case of n expenditures and has proved the following theorem: 





Theorem. If an Engel Curve is given by x; = a; t? + b.t+c,, 4 
t= 1, ..., nm, where t is the total income and x. is the amount spent - 
spent on a particular item then the utility function, U, is of the fom < 
Us K(- 2x? - x,*) - Ax,, ..-, %,) + F where K is a arbitrary positive 5 
constant, F is any constant, and Ax,, vee, XD 2S arbitrary positive e: 
definite quadratic form in (x,, aes x). 

4. Summary. As a general mle the expenditure on an item or on a group 1 
of closely related items is a quadratic function of the total expenditure : 
or income. 0 

The situation for two variables, savings and al! other expenditures, 
is discussed in detail. The Engel Gurve is a parabola. The utility s 
function, if quadratic, is determined except for an additive constant E 
and a constant multiplier, by this parabolic curve. The indifference b 
curves are arcs of ellipses and the quadratic utility function 1s convex. 0 

The n dimensional! case permits additional elements of arbitrariness S 
in the utility function. n 
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SYSTEMS OF EQUATIONS, MATRICES 
AND DETERMINANTS 
(Conc luded) 
Olga Taussky and John Todd 
CHAPTER II 


The numerical treatment of problems in this field is not entirely 
straightforward, and requires particular care even when we are handling 
problems 1n quite a moderate number of unknowns e.g., for about 10. 
This whole subject 1s at present under active investigation in view 
of the development of high speed artomatic digital computing machines, 
by means of which 1t 1s possible to contemplate the solution of problems 
in which the number of unknowns is of the order of hundreds. These 
investigations are in various directions. Me of these is a re-examin- 
ation of old methods and includes their history, classification and 
unification, together with detailed studies of special methods and 
attempts at the eradication of the superstitions and the justification 
of the suspicions which are still mfe in the rather primitive field 
of numerical analysis. A second is the devising and examination of 
newer methods, such as the gradient or finite iteration methods. A 
third direction 1s the study of methods appropriate for special systems, 
such as those which arise by the discretization of differential equa- 
tions. 

We shall show some of the advantages and disadvantages of the class- 
ical methods, indicate some of the newer methods, and show some of the 
difficulties which can arise. This will be done mainly by discussions 
of special numerical examples. 

For up -to-date accounts of this aspect of the subject, reference 
should be made to Proceedings of a Symposium on Simul taneous Linear 
Equations and the Determination of Eigenvalues’. The article in this 
by G. E. Forsythe contains a comprehensive bibliography of the first 
of the two topics. For a more detailed account of methods for the 
solution of equations and the inverting of matrices, illustrated with 
numerical examples and directed more towards those who do not have 
access to high-speed equipment, we refer to the report by L. Fox*. 


Il. 1. SOLUTION OF SYSTEMS OF SIMULTANEOUS EQUATIONS 
Methods of solution can be classified as direct, or as indirect 


(or iterative). In the first, in theory, we obtain the exact solution 


'To be published as National Bureau of Standards, Applied Mathematics 
Series, Vol.29, 1952. 


2To appear in National Bureau of Standards Journal of Research, 1952. 
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after a finite number of steps. In the second we produce an infinite 
sequence of numbers which converge to the exact solution. In exceptional 
circumstances the sequence may be stationary, i.e. all terms are equal 
after a certain stage. There have recently been developed methods 
which, in theory, are stationary in general. In practice however, 
because of the fact that our calculations have to be approximate, they 
are not stationary. 


DIRECT METHODS 


l. a. Elimination Method 
This is one of the oldest methods, and, it will be seen, one of 
the best, at least for general systems. Consider 


(1) R, = 12x - 3y + 22 - % = 0 Check sum = -85 
(2) R, = -3x - By + z - 68 = 0 Check sum = -78 
(3) Ry = x + dy + 6z2- 3=0 Check sum = 6 


Eliminate x by subtracting 12 times the third equation from the first 
and by adding 3 times the third equation to the second. We get 


(4) -27y - 70z - 60 = 0 Check sum = -157 
(5) <-2y +19z-77 = 0 Check sum - 60 


" 


Eliminate y by multiplying (5) by 27/2 and subtracting from (4). We 
find 


653 1959 
(6) -—-z+-——=0 Check sum = +653 
2 2 
which gives z = 3. We now substitute this value of z in (4) (or (5)) 
to get y = -10 and then put these values in (1) (or (2) or (3)) to 
get x = 5, 


In this, as in all numerical work, a description of a proposed 
solution is not complete unless some sort of checking system is in- 
corporated. This should reveal errors as soon as possible after they 
occur and not at the last stage, if at all. There is available in this 
case, the following method which can be applied with minor modi fica- 
tions in most of the manipulations in this field. We shall not mention 
it explicitly again. 

We carry an additional column in our wrk sheet as the sum of the 
numbers in that row. After writing dow our system we compute thie 
numbers -85, -78, 6 labelled “Check sum” above. We perform on these 
sums the same operations as we do on the equations. Thus we take 
-85 - 12x 6 = -157 and -78 + 3x6 = -60 ad compare these with the 
sums of the coefficients in (4) ad (5) which are -27 - 70 - 60 = -157 


and -2 + 19 - 77 = -60. Proceeding, as a check sum we compute 
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-157 - [x( -A0)] = +653 which agrees with the doefficient sum in (6): 


653. 1959 


—s > = +653. A final check might be the substitution of the 


~ 


values of x, y, z in the two equations which were not used in the 
determination of x. 


1. 6. Determinantal Solution 
This 1s often called Cramer’ s Rule. Applied to our example it gives 








x -y z =} 
ee eae Chg 
}-3 2-9} =| 12 2 -% | =| 12 -3 -%| =|12 -3 2 
}-8 1 -68 | -3 1 -48 | -3 -8 -68 -3 -8 1 

26 - 3] 16 -3 1 2-3 11 2 6 
1.€ 

x ~y z 1 

3265 : 6530 : 1959 "653 


mving x = 5, y = -10, z= 3. Thus the problem is reduced to the evalua- 
tion of determinants. 

The evaluation of a determinant of order n from its explicit defini- 
tion becomes rapidly more tedious as n increases. For there are n! 
terms in the expansion and each involves n factors; we have therefore 
to carry out about n! xn nultiplications as well as about n! additions. 
A more efficient method of evaluating determnants is therefore essen- 
tial. We shall show how it is possible to transform a determinant by 
means of the transformations of the form discussed in [. 2, which do 
not alter its value, into one whose expansion contaliis but One non-zero 
term. Take the determnant 


12 -3 2 
wei de 8 ht 
1 y 2 6| 


Denote its rows by ry> To, T,- Subtracting twelve times r, from r, ad 


adding three times r, to r, we obtain: 
| 0 -27 -70 | 
D = | 0 -2 19 
io 2 6 





by 27/2 and subtracting from r, we obtain 


Multiplying r , 


2 
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-653 
2 
0 -2 19 


1 2 6 








The expansion of this determinant contains the single non-zero term 


653 
-(1) (-2) (-—) = -653. 


What is the expense of this operation in the general case? Again 
neglecting additions we see that the first stage contains one division 
and nX(n- 1) multiplications. Neglecting the divisions we see that 


we have in all about > n(n - 1) = _* multiplications. This is a 


considerable improvement on the n!x n multiplications, even for smal] 
values of n, 

This method, which we see is essentially equivalent to the elimination 
method, is recommended for the evaluation of determinants of general 
type. It also seems reasonable to discard the determinantal method 
of solution, even if the determinants are evaluated by the method just 
described, in favor of the elimination method. This method, as applied 
to determinants often goes by the name of Chidé’s method of pivotal 
condensation. The pivots in our example are the numbers 1, -2. There 
is, of course, plenty of freedom in the choice of the pivots (i.e., the 
order in which we eliminate the variables). The most efficient choice 
is being investigated. 


l. c. Finite Iteration Scheme 
The standard method for fanding the center of an ellipse 


(1) ax? + Qhxy + by® + Qgx + 2fy + c = 0, ab > h? 
is to solve the system of equations 
(2) ax + hy + g=0, hx + by + f = 0. 


Its solution gives the center of (1) for any c. We shall now reverse 
the process: given a system of equations of the form (2) we shall 
construct the center of the famly of ellipses (1). We use the fact 
that the chord of contact of a pair of parallel tangents passes throug) 
the center. Let 9, be any direction. Denote by |, the chord of contact 
of the tangents in the direction 9,. We denote by 9, the direction 
of 1, ad by 1, the corresponding chord of contact. The intersection 
of the lines 1), 1, is the center. We reach it as follows: Let S, be 
any point. Proceed from S, in the direction 9, until we meet l,, at 
S, say. Proceed from S, in the direction 9, until we meet 1, at 5S): 
Then S, is the center. 
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The method can be extended to higher dimensions; for instance, con- 
sider the enveloping cylinders, in three directions, of a ellipsoid. 
The three planes of contact intersect at the center. This can now 
be approached in three steps. 














Fig. 1 


We shall, for simplicity, in the diagram work out a two-dimensional 
example. Compare Fig. 1. To solve 


( 3) x-y2O0,-x+ 2-120, 


we determine the center of the ellipse 


x* - Qxy + 2y? - dy = 0, 

We take S, = (0,0) and observe that the residuals in (3), i.e. the 
values of the left hand sides of the two equations for x = 0, y = 0, 
are 0, -1. We choose the direction of the residual vector (0,-1), 
in this case the vertical direction, to be 9, The chord of contact 
1, is indicated in the diagram. We find the point S, = (0,%). Then 
we obtain the chord of contact 1, of tangents parallel to 1,. We then 
proceed from S, along 1, until we meet 1, at S, = (1,1), the center 
required. 

The general algebraic description of this method for solving Ax = 6 
is the following. We choose any S, and observe the residual vector 
R = AS, - 6 (if R, = 0 we are finished). We choose Z, = R, and define 
successively 





MATHEMATICS MAGAZINE (Nov. -Dec, 


Zin, = Ry - 5AZ, + 0,2;, Ry,, = Ry - a, AZ; 
where the a;, 6; are certain scalars: a; is chosen to make S;,, the 
proper distance from S; in the direction Z, and then 6; is chosen to 
make Z;,, parallel to the planes of contact of all preceding Zi. 
We have to take 
2 . tt 2 

a, = R,/ZAZ,, 6, = Ry, /R;. 
It 1s found that R;,, 1s indeed the next residual and that R, is 
certainly zero in theory. In practice R. will be small, and it may 
be necessary to repeat the process. 


INDIRECT OR ITERATIVE METHODS 


1. d. Relaxation 


We now retum to the 3-dimensional example of l.a. The solution of the 
equations (1), (2), (3) 1s accomplished by choosing values of x, y, 2 
to make the residuals R,, R,, R, zero theoretically or very smal! 
practically. We consider the effect of unit changes on the values of 
the residuals and obtain an operations table 


R R R 


1 2 3 
12 -3 
-3 -8 
2 l 


The matrix here is the transpose of the matrix of the system. The relaxer 
tion process, in its most naive for, starts off with arbitrary values 
of x, y, z and at each stage liquidates, as nearly as possible, by 
altering one variable, the largest residual. Considerable virtuosity 
in the art can be achieved by practice, e.g. instead of altering a 
single variable several can be altered in a “block-relaxation”’. This 
may be suggested internally from the behavior of the residuals or 
extemally from some knowledge of the synmetries in the physical problem 
which gives rise to the numerical problem. 


R, R, 

-96 -68 
0 -92 
33-4 

-3 
3 8 -2 
0 


Check 0 
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Starting with a guess 0, 0, 0 we obtain the residuals -96,-68, -3. 
Looking at the operations table we see that the most economical way 
to liquidate R, = -96 is to change x by 8 thereby al tering R,, R,, R, 
by 96, -24, 8. We obtain the second row above. Now R, is the largest 
residual and we can liquidate it approximately by a change of -]l in y, 
which means altering R,, R,, R, by 33, 88, -22. We thus get the third 
row above. We now reduce R, by a change of -3 1n x, the R, by a change 
of 3 in z and then R, by a change of 1 in y. At this stage, in view 
of our specially simple example, we get zero residuals. We check our 
work by adding up the changes in x, y, z to get x = 5, y = -10, z= 
and check that the residuals are acwally 0, 0, 0. 

In practical cases we rarely reach an exact solution. What usually 
happens 1s that the residuals are reduced by an order of magnitude. 
If this is not sufficient we change the scale in our residuals and the 
variables by a factor 10 or 100 ad begin again from this first approx- 
imation. This can be repeated until a satisfactory solution is obtained. 

Thus described, relaxation is a paper and pencil method and we need 
only work with small integers. The process can be mechanized, but it 
would be difficult to codify all the tricks of the trade. We point 
out that the method is particularly effective when the matrix has a 
dominant main diagonal. 


l. e. “Seidel” Iteration Scheme 

We again take the three-dimensional case. We begin with any guess 
at the solution, say . * l, y= a 8, * 1. We obtain first a revised 
value for x by substituting the old values of y and z in the first 
equation and solving for x. This gives x = 8.08. We then use this 
value of x and the old value of z in the second equation to obtain 
a revised value of y:y, = -11.40. We then use x and y in the third 
equation to get a revised value of z = z, = 2.95. We improve this 
first approximation 8.08, -11.40, 2.95 in the same way. And so on. 
A few stages of the process are indicated below: 


l 12x - 3+2 = 9% : 8.08 
l -24.24 - By + 1 = 68 : -11. 40); 
8.08 - 22.80 + 6z = 3: 2.95 


12x + 34.20 + 5.90 = 96: 4.66 5.0 
-13.98 - By + 2.95 = 68 -9.88/, ...| 10.01 
4.66 - 19.76 + 62 = 3 3.02!) 3.00 


This process is convenient both manually and mechanically. Convergence 
is assured if the matrix is positive definite. 


1. f. Gradient Methods 


We return to the example discussed in l.c. We observe that the 
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solution to the (consistent) system 
x- y=0, -x + 2y = 1 
is the point which minimizes the quadratic form 
€(x,y) = (x - y)? + (-x + 2y - 1)? 


the minimum being zero. 








Fig. 2 


We take an arbitrary point e.g. (x,,y,) = (0,0) as our initial approx- 
imation. Compare Fig. 2. We wish to proceed as rapidly as possible 
to the point where €(x,y) is minimum. Now €(x,y) decreases most rapidly 
at (x,,y,) in the direction of the inward normal to that ellipse of 


the family 
(1) €(x,y) = constant 


passing through the point. We therefore begin in this direction; it 
seems reasonable to continue moving in it until we have reached the 
minamum value of €(x,y) on this nomal: this will occur when the nomal 
becomes:a tangent to a member of the family (1). With the choice 


( de 
(x ¥5) = (0,0) the direction is ail 


Ox 
0,0 


($-| | (-2,4). Ou 
oe » eee 
Wi, 4 


next approximation is 
me By Ste Fy * Hy + 
where r is such that 
(x, - y,)* + (-z, + 24y, - 1)? - k 


is a perfect square (in r) for a suitable value of k. We find 
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k 9 =3 10 ( ’ 9 
c—- fF = —, s-, El x P = - 
34 if *s" Se 7 "34 


We then proceed from (x,,y,) in the new direction of steepest descent: 
A 63 
a = (2,1). 
17 17 
25 25 2)" 
We find (x,,y¥,) = ==. =F e(x,,y,) = =a . In a similar way 


we fin d 


' y = {805 a) ' 
x ’ y = ental aver we E(x ’ » 
= (1156 1156 3°73 


and 


sang (RE, ES. a | 


, {, ECX 
1158 1158) (34) 


The approach of (x,y,) to (1,1) is indicated geometrically in the 
diagram. The points GS, S,, o- lie alternately on the lines y = x 
and 8x - l3y + 5= 0. 

The question as to whether it would not be advisable to go beyond 
or stop short of the point of tangency has been investigated in special 
cases, but no general decision has been made. 


II. 2. INVERSION OF MATRICES 


The solution of the system of equations Ax = 6 can be obtained by 

application of the matrix A'' to the vector b, an operation involving 

multiplications. If it 1s desired, as often happens in practice, 

solve the system Ax = 6 for various values of 6 (A being undiamyed), 
then it is frequently advantageous to invert the matrix A, once for all. 

Many methods are available for this inversion. They can be classi fied 
as direct and indirect or iterative as in the previous section. Many 
of the methods for the solution of systems of equations can be adapted 
to the more general problem of inversion. 


2. a. Cholesky’s Method 

This is essentially based on the elimination method. In its sinplest 
form it applies to a symmetric matrix. We shall discuss the inversion 
of a matrix A by using its representation as the product of an upper 
triangular matrix U and its transpose. Thus | 


12 3 fa 0 a b a® ab ac 
2 3 44=u'U=|b6 d O10 d ab b*+d* be+de 
3 4 4) c e filo 0 ac be+de cz+e%s f*) 


With this assumption we find in turn a? = 


\ 
4, 
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ac = 3, c= 3; b* + d*= 3, d= i, be + de= 4, e = 2i and c* + ec? + f* 
= 4, f = t, where we have made arbitrary decisions on the ambiguous 
signs of c, d and f, which are determined by their squares. It is to 


be observed that no simultaneous equations have to be solved. 
We next observe that 


A’ = (u'uy”' = (u')yqu’)' = uo (U')! 


wt at 
and that U- is easy to find. In fact U is also an upper triangular 
matrix which can be computed from the relation 


2 3 


om 


¥2 


1 
0 
0 


ap 
l 0 0 € 
0 ijJOW O 
from which we obtain successively a 
€ = 21, DP= -t. Yence we have 


1 2i -vyrl 0 0 -] 
A'=]/0 -i 2i]/;/2i -i 0O 
0 O -i}l-i Qi -i -l1 2-1 


The method is sometimes called, for obvious reasons, the square 
root method. We observe that in it we may introduce surds and complex 
numbers but that these do not appear in the final result - inversion 
is a strictly rational process. 

The second remark is this, that if we apply the Cholesky process 
to the solution of equations we can reduce our wrk somewhat. We replace 
the equation Ax = b by U'Ux = 6b and this we can replace by y = Ux 
and U'y = b. We therefore solve first for y and then for x. In the 
case when A is the matrix under consideration and b the vector (14, 
20, 23) we can arrange the results of the work compactly as follows: 


l 3 14 ] 14 1 
2 4 20 81 2 
4 23 3u 3 

b j y x 


We also note, that, in addhtion to linear checks of the type mentioned 
in 1.2, we can use a quadratic one. We have 


x'b = x'Ax = x'U'Ux = (Ux)'(Ux) = y'y 
In the present example we have 


14x1 + 20x2 + 23x 3 = 123 and 147 + (Gi)? + (31)? = 123. 


2.b. Iteration 


It is well-known that the sequence of real numbers defined by the 
recurrence rel ation 
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Xue, * %,(2- Nx.) 


converges to N' for suitable x,, for any real number N #0. It can 
be shown that, for suitable X,, the sequence of matrices X, defined by 


oe = X (21 - AX,) 
converges to A'', it being assumed that A is non-singular. As an example 
take the matrix A to be that of 2.a. If we take 


-3.0 4.1 -0.9 
4.1 5.3 1.9 
-0.9 1.9 -1.1 


we obtain 


-4. 26 4.14 -0.6 9976 3.9864 -1.0138 
4.14 -5.06 1.94 9864 -4.9896 2.0104 
-0.86 1.94 -1.06 -1.0136 2.0104 -0.9896 


This process 1s useful for the improvement of approximate inverses. 


Il. 3. CHARACTERISTIC ROOTS OF MATRICES 


It has already been remarked that the problem of determining tlie 
characteristic roots of a symmetric matrix is essentially the problem 
of finding the principal axes of a quadric surface. These pure mathe 
matical problems underlie many physical and engineering problems concem- 
ed, for instance, with the theory of vibrations. 

There are two different problems in this section; according as to 
whether we wish to determine a few of the dominant characteristic 
roots (i.e. those with largest moduli), or whether we wsh to find all 
the roots. The first of these problems is the simpler. The following 
comments refer to the second problem. In the first place, our remarks 
about the evaluation of determina@ts indicate that a direct expansion 
of the determinant, and the assembly of terms into the characteristic 
polynomial 1s impractical even if the solution of a polynomial equation 
could be obtained readily (which is far from being the case). Again, 
let us assume available a practical method for obtaining a dominant 
root and a method for obtaining a matrix (of one lower order) having 
the remaining roots of the original as its characteristic roots. It 1s 
tempting to say that repeated application of this process will produce 
all the characteristic roots. In point of fact, however, the dominant 
root will only be obtained approximately, and the new matrix will 
therefore be inaccurate on two counts: due to the inaccuracy of the 
root and due to the inaccuracies inherent in the numerical calculation 
involving this approximation to the root. Very soon indeed all sig- 
nificance will be lost. Among the more satisfactory solutions to the 
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problem is one due to Jacobi: it is described in [I. 3. c. 
We shall confine our attention to the case of real characteristic 
roots. This covers the case of symmetric matrices completely. 


DOMINANT CHARACTERISTIC VALUES 

3.a. Iteration 

Choose an arbitrary vector v normalized in some way, e.g. to have 
one of its coordinates unity or to have the sums of the squares of 
the coordinates unity. Apply the matrix A repeatedly to the vector 
yl), expressing each product vector as a scalar multiple of a vector 
in the chosen normalization. Specifically, if v't) is normalized then 
we define u''*!') by the equation 


(o) 


Av'*) = ylttt) Ctr 


where v''*') is normalized. It can be shown that if A has a single 
dominant root then these multipliers tend to this value ad the (nomal- 
ized) vectors tend to the corresponding (normalized) characteristic 
vector. 


Let us consider the case 


0.2 0.9 1. 32 
= i-11.2 22.28 -10.72 
-§.8 9.45 -1.94) 


and choose v'° 


) s (1,0,0) and normalize by making the first coordinate 
unity for simplicity. We obtain the following results: 


Av'°) = (0.2,-11.2,-518) = u!')v'"') where u'') = 0.2, 
vil) = (1, -56, -29) 

Av'') = (-88. 48,-948,-478.74) = u'?)v'?) where u'?) = -88. 48, 
vi?) = (1, 10.7143, 5.4107) 

Av'?) = (16.9850, 169. 5119,84.9534) = u'*)v'>) where u'?) = 16.9850, 
v' 3) = (1, 9.9801, 5.0017) 

Av'3) = (15.7843, 157.5384,78.8086) = u'4)v'4) where u' 4) = 15.7843, 
v4) = (1, 9.9807, 4.9928) 


Av‘ 4) = (15.7731, 157.647 2,78.8316) = u'5)v'%) where u'®) = 15.7731, 
v'5) = (1, 9.9947, 4.9979) 


Av'5) = (15.7925,157.9044, 78.9540) = u'®)v'6) where u'®) = 15.7925, 
v'6) = (1, 9.9987, 4.9995). 
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It can be verified that the exact resuJts are A, = 15.8 and 
v, = (1, 10, 5). 

The justia fication of this process is simple. It 1s know that, common- 
ly, a matrix A has n different characteristic roots A; adn distinct 
characteristic vectors c; which are linearly independent and which 
therefore span the whole space. An arbitrary vector y!?) 


in the form 


cen be expressed 


io? « 
Vv = 2a, c;. 


2, «++, nm we have 


yim) = gMylo) 2a Ale; 


and from this, if A, >A, > *** >A, we have, for sufficiently large 
n (depending on the separation of the A’s) that 


yim) 2 grylo) ani cy 


From this the statements made above follow: that yin) is approximately 
a multaple of c, and that the ratio of corresponding components of 
v'™) and v'"~') is approximately r,- 

3.b. Relaxation using Rayleigh’s Approximation 

For this method we recall the fact that the dominant characteristic 
root A, of a symmetric matrix A is the upper bound of the Rayleigh 
quotient 


xAx' 


‘ 
xx 


R(x) = 


taken over all vectors x and that this bound is assumed wien <x is the 
corresponding characteristic vector, c,. It is also know that the 
difference between R(x) and A; is comparable with the square of the 
distance between x and c,, when x is near enough to c,. We apply the 
Relaxation method to systems of equations which are not quite com- 
patible since they involve an approximation to the characteristic 
root sought. 
Take the case when 


3.5000 0.750 1.299 
= 10.750 1.625 1.083 
1.299 1.083 2.875 


We guess v'°) = (1,1,1) and estimate the dominant characteristic value 
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3.500 + 1.625 +2.875+ 2[0.750+ 1.299 + 1.083) 
ul!) = R(v) s-_— —s 
12 + 12 + 12 





14. 26 
= : = 4.755 





We consider the homogeneous system of equations wth matrix A-u"')J ie 


-1.255 750 §=1.299 
.750 -3.130 1.083 
1.299 1.083 -1.880 


We relax according to the following scheme 


se i yet see % Residuals .794 -1.297 . 502 
-0.5 . 419 .268 -.040 
-0.1 . 289 . 160 .148 


Sum: x=] y=0.5 z= 0.9 Check: .2891 . 1597 -1485 








We now use v''?) = (1, 0.5, 0.9) as a revised approximation to 


the characteristic vector and re-estimate the characteristic values 
ul?) s RO v'')) i 





3. 500x1+1.625x0 . 25+2.875x0.81+2(0.750x1x0. 5+]. 299x] x0.9+1. 08 3x0. 5x0. 9] 
re 1+ 0.25+0.81 





= 4.999 


We now consider the revised system of equations with matmx A - ul?) 7, 
1.e 
-1.499 .750 1.299 
0.750 -3.374 1.083 
1.299 1.063 2.126 








We relax, beginning with (x,y,z) = v!") thus: 
zs*<j] y 26.5 2296.9 Residuals .045 .038 -.071 
-0.3 .006 .006 -.007 
-0.003 - 002 .003 -.001 
Sum: x= 1 y= 0.5 z= 0.867 Check: .0020 .0020 -.0010 


Thus we can say that 4.999 is @ approximation to the dominant character- 
istic value of the matrix and that the corresponding characteristic 


vector is (1, 0.5, 0.867). 








to 





ES 
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ALL CHARACTERISTIC ROCTS 


3.c. Jacobi’s Method, 

The method to be described is based on the fact that by an orthogonal 
transformation from variables x,y to variables x‘, y' which we can 
describe in the form 


[* - cos 9, sin 9 ‘| 
y' -sin 9, cos 9)ly 
we can express the quadratic form 


ax* + Qhxy + by? 


as a sum of squares 
Zz #2 
Ax + By‘; 


we have only to choose 9 such that tan 2°=2h(a- 6) '. 
We take the matmx 
2.879 -0.841 -0.148 0.506 
-0.841 3.369 -0.111 0.380 
-0.148 -0.111 1.216 -0.740 
0.506 0.380 -0.740 3. 536 


We observe that the largest off-diagonal element is -0.841. In view 
of the result just quoted it follows that if 


tan 26 = -2 x 0.841/(2.879 - 3.369) and 


cos 9 sin? 0 O 
-sin 9 cos 9 0 
Ty = 
0 0 0 
l 


oOo - © 


0 0 


then the element in 7,AT, corresponding to -0.841 will be zero, or 
rather, very small in view of the fact that our calculations are not 
made exactly. We find, in fact, that 

2.250 0.001 -0.186 0.43 

0.001 4.000 0.002 -0.002 

-0.186 0.002 1.214 -0.740 

0.636 -0.002 -0.740 3.534 


A. = T, AT, 


The largest off-diagonal element is -0.740 and we reduce it by trans- 
formation by Ty where 
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1 O 0 )) 

aa 0 0 
Ty = 

0 O cos ? sin 

0 Q -sin DP cos 2 


and where tan 20 = 2x(-9.74!)/(1.2146 - 3.536). We find 


2.250 0.001 0.000 0.66 

0.001 4.000 0.001 -0.002 
0.000 0.001 1.004 0.001 
0.663 -0.002 0.001 3.755 


A, = TjA T= 


We next reduce the element 0.443 by transformation by Ty where 


cos¥Y 0 0 sin W 

0 1 O 0 

N° 0 oO1 0 
l-siny 0 0 cosy 


where tan 2/ = 2x0.6463/(2.250 - 3.755). We fin} 


4.905 -9.002 0.000 -0.001 
-0.002 4.000 0.001 -0.002 
0.000 0.001 1.004 0.0900 
-0.001 -0.002 0.000 1.999 


Att 


suggesting that the characteristic roots, wiich are invariant under 
the above transformation since tle 7’ s are orthogonal, are aproxmatel y 
4.005, 4.000, 1.004, 1.999. 

The proofs that this process is convergent and that it is a practical 
one are not hffieult. 


II.4. ROUNDING-OFF ERRORS 


The question of the effect of rounding errors in calcul ations is 
rather difficult but of extreme importance in practice. We shall con- 
fine our attention to the case of matrix inversion. What is required 
in this case are bounds for the difference (in some sense) between the 
true inverse of a matrix and the result obtained by carrying out a 
particular process, working always to a fixed number of decimal (or 
binary) places. These bounds should be readily obtainable and may, 
of course, be estimates, not necessarily very precise. The bounds 
may be absolute ones, or probabilistic ones: in the sense that n€ is 

n 


the absolute error in the sum 2 a, due to errors Ei le, | < €, in 
t=! 





a de cl 
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a, while the probabilistic error 1s estimated as 6&/n/12. 


No satisfactory solutions to these problems exist and the results 
now available are too complicated to diseuss here. To underline the 
importance of these problems we shall discuss a particular numerical 
example. 

Consider the application of some of the methods just described to 
the case of the solution of the system of equations 


10x + Ty + &z+ Tw 32 
Tx + Sy + 62+ Sw = 23 


8x + fy + 10z + Aw = 33 





7x + Sy + 92+ 10w = 31 


which has the solution x 
by T. S. Wilson. 

If we eliminate first y and then w and wrk only to one place of 
decimals we obtain 


" 
< 
" 
N 
" 
e& 


»= 1. This example was constructed 


O.lx + 0.3z = 0.4 


0.4x + 1.22 


1.6 


These equations are not independent so that elimination of x also 
eliminates Zz, 
If we carry three decimals we find 


0.060 z = 0.04) 


which gives a poor determination of z = 1.07? 
If we apply the relaxation process, we can arrive at the following 
results, where R,, R,, R,, R, denote the successive residuals 
1 R, R, R, 
16 .6 “Jem “2.9 3.1 0.1 -0.1 -0.1 0.1 
2. 34 0.18 0.45 1.21 0.01 -0.01 -0.01 0.01 


1.136 0.918 0.965 1.021 0.001 -0.001 -0.001 0.001 


x vy z w R 


If we apply the Seidel process we obtain tle followng sequence of 
approximations: 


0) (3.20) (2.44) (1.98 (1.70) (1.5 
0} [0.12] [0.18 |o.21] |0.22 oe 
o|* Jo.e7]’ [1.061 ° }1.28)° 1.30] * Ja. al’ 
0) lo.20) (o.35) lo.46) lo.ss) Loven 
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The behavior of this system is blamed on its lack of “condition”. 
Various attempts have been made to measure this lack of condition. 
The most primitive of these is the smallness of its determinant (in 
comparison, for example, with the individual terms in its expansion). 
More satisfactory measures have been introduced. Among these are the 
M-, N- and P- condition numbers defined by 


MA) = n R(A)R(A') where ® = (C) denotes max le; ; | 
t,J 


MA) = n°' norm A nom 4’! 


P(A) = A/u where A, & are the greatest and least of the absolute 
values of the characteristic roots of A. 


When these numbers are large, trouble is to be expected. The actual 
values of the numbers in the case under discussion are respectively 


27 20, Tale 2984 
wiich are to be compared with the values 
5.60, me 22, 9.47 


for the matrix 


which 1s typical of those which arise by the discretization of dh ffer- 
eitial equations. 


Estimates for the error in matrix inversion can be given in terms 
of the condition numbers but it is clear that these numbers cannot 
be calculated wmthout a knowledye of the inverse (in the first two 
cases) or of the characteristic roots (in the third). 


ERRATA: 


Page 13, Sept.-Oct., 1952. The equaticn 
l 


oe L+k 

4 * Tay ((-1)** A.) 
shou] d have read 

A°' s 


1 
TAT ((-1)**A,,)'. 
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TEACHING OF MATHEMATICS 
Edited by 


Joseph Seiclin and C. N. Shuster 


This department is devoted to the teaching of mathematics. Thus artic'es on 
methodology, exposition, curriculum, tests anc measurements, and any other 
topic related to teaching, are invited. Papers on any subject in which you, as 
a teacher, are interested, or questions which you would like others to disciss, 
should te sent to Jeseph Seidlin, Alfred University, Alfred, New York. 


IDEAL PREPARATION IN MATHEMATICS FOR COLLEGE ENTRANCE* 


J. Seidlain 


The tatle 1s not of my making. But for students and teachers of 
Nathematacs 1t probably conveys the intended meaning. In the study of 
ideals or limits we are generally concerned with the behavior of the 
variable or variables approachiny our ideal or limit. Specifically, 
then, our topic, sheared of symbolism, beconwes: “How may grade teachers 
an! high school teachers do a progressively better job of preparing 
their pupils for optinum learning in courses in college mathematics? 
That’s it; nothing more and, perhaps, nothing less. 

On April second T wrote to forty of my friends and colleagues, as 
follows: 

‘Dear Professor X: On May 3rd I am to telj an assembled multitude 

of secondary school teachers of mathematics something about “Ideal 

Preparation in Mathematics for College Entrance”. It occurs to me 

that I should serve “the cause” better if, in addition to my ow 

judgment, I would gather the judgments of the “workers in the fheld”. 

Would you write me as much or as little as you please, on what you 

cousider the ideal preparation in mathematics for college entrance? 

May I quote any part or all of what you say?’ 

I got a wonderful response and I could do wrse than just reproduce 
what they said. 


“relevant di- 


At this point J wait to digress a bit, but itis a 
zression’’. 

On December 19, 1931, I read a paper before the Mathematics Section 
of the New York Society of Experimental Research in Education. A year 
later, December 1932, this paper was reproduced in the Mathematics 
Teacher under the nasty title of “ The Contribution of High Schoo! 


Teaching to Ineffective College Teaching”.** Then, as now, I had sought 
*Paper read at the Seconc Annual Meeting of the Association of Mathe- 


matics Teachers of New York State, May 3, 1952, at Syracuse, N.Y. 
**The Mathematics Teacher, December 1932. (Vol. XXV, Number 8) 
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and received the help of teachers of college mathematics. Twenty-one 
years ago T grouped all the criticisms of high school teaching and 
teachers under twelve headings. The following is a fair sampling: 

1. “They (the students) come to us jaded, worn out, weary of math- 
ematics” , 

4. “They do not know the principles underlying the simplest opera- 
tions.” The greater rumber of our students become involved in the so- 
called advanced work because 

e. They are generally weak in the application of the laws of 
exponents (literal, fractional, negative). 

f. They are extravagant in the use of the terms ‘‘transpose’’ 
and “ cancel ’”’. 

6. “Not one in ten is conscious of the positional schere of our 
nun’ er system. That is why they experience such difficulty in vorkiny 
with detached coefficients” . 

1]. ‘They are happiest when a ‘proof’ (demoistration) is omitted”. 

Primarily, the high school teacher is not blamed for failing to 
impart a sufficiently large number of facts to his pupils, but for 
apparently failing to inspire, or keep alive in, the pupil the nore 
genuine values presumably obtainal le through the study of mathematics. 
The implicit query raised in most instances is: “Of what conceivable 
value are more or Jess perfected schemes for imparting facts - facts 
soo: forgotten - if both the founcation amd the superstructure of these 
facts are so fatally impaired?” 

What new or different criticisms do we find after a lapse of twenty- 
one years? Before I answer this question specifically, I wmt to make 
a rather sisnificant and encouraging generalization. Today, as contrasted 
with twenty-one years ago, the college teacher seems to be more keenly 
aware and appreciative of the high school teacher’s enormous burden 
of trying to meet all the needs cf all the children through mrmeular, 
co-curricular, and extra-curricular tasks. Thus, while criticisms aimed 
at the latest product of high school teaching are as plentiful as ever, 
I detect an attituce of chamty, a note of sympathy, and even a modicum 
of humility. 

And now let us get a glimpse at the 1952 edition of “helpful hints’ 
to the producers of freshmen in college mathematics. 

1. (amd most prominent) Fundaverntals and understanding. 

(a) ‘'. . . IT wish the young people coming to college really under- 
stood the fundamental ideas of algebra and trigonometry 7. 

(b) “. . .the ideally prepared student should understand. . 
whatever topics he may have studied in high school . 

(c)“. . . in the preparation of students more stress should be 
laid on understanding .. .”’ 

(d) “. . . unless a student’s mathematics has a strong foundation, 
any super-structure you try to build. . . is a pretty flimsy 
affair... .”’ 
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(e) “. . . the number of topics covered is not as important as 

a thorough understending of U.ose that have been covered...” 

In all but two of the replies “understand” or “understanting” is the 
key word, 


2. Subjects. 


Arithmetic an? Elementary Algebra lead the list. Plan» Ceometry is a 
close third, All of the other secondary school branc'.es of mati.ematics 
are mentioned at least once and for our lower and upper limits, respec- 
tavely, we have the following: “... T think it as relatively unimportant 
what subject matter is chosen”, and“... this world imply that the high 
schools have given their college matriculants competence in al! branches 
of mathematics tlhro:gh the calculus”. 

However, pretty generally, and stated more or less kindly, we find 
the plea that analytic geometry and calevlus shovld be judiciously 
imored in tlhe high schoc} course of stu:ly. Under a less frequently 
mentioned “hands off” policy we fird even advanced algebra and sclid 
geometry, though in the case of these two, recognition is given to the 
fact that some engireering sclhcols require one or bot! of them. 

All ain al], at is not difficult to summarize the domnant position of 
college professors of mathenatics wth respect to the prerequis.tes for 
college courses. Tt 1s not how much mathematics is covered, but how 
fundamentally wel! at 1s learned and retained by the pup:!s. 


3. Boners. 


There seems to be no Jimit to the ingenuity displayed by students in 
juggling figures or msapplying teclniques or perverting principles to 
suit. their own dim and veiled ends. No teachin, or teacher, however 
poor, should he held responsille for this seemingly countless variety 
of pupils’ boners. Stall, most of us feel that better teaching would 
reduce the strain on our sense of humor. 

Vany of these “intangible errors’’ are familiar to all of us. Yet 
every now ac then f min across a new one. Why doesn’t someone collect 
them and publish them as a sort of chamber of horrors? Several of my 
respondents presented me with a few samples, viz. 


l ] ] \ Proof of a special kind 
4 ‘ 4 | of cancellation 








f 
x - af f(x) Proof of the RemainZer 


f(x) - fla) Theoren 
+ fla) 
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1 
a 


and many others. 


2 


1 
b a + b 


+ 





4. Quotable quotes. 


“ 


-+» I feel that the bigh schools in New York Stare 4b a fairly good 
job of preparation in the field of Mathematics’’. 


rf 


+. every attempt should be made to make the student fec! that the 
material studied in one course can be transferred to another course”. 


“A student who knows the definition of the quotient of two numbers wil! 
not wonder why he can’t divide by zero”. 


“Many college freslimen cannot read even a daily newspaper”. 


“The most frequent and important practica! use of trigcnometry today is 
not in the solutions of triangles but in the study of wave motion and 
other periodic phenomena”. 


“In view of the appalling state of mathematica] ignorance in which 
students enter and too often leave our colleges - 


“The best definition of a circle that a class of thirty fres!men could 
produce was ‘a circle is when something goes around’ ”’. 


“Whether the basic difficulty is with the teachers or with the syl! abus 
is something that perhaps could be thrashed out at the Syracuse meeting”. 


... they (the students) are so much disturled and frustrated by the 
arithmetic and elementary algebra entailed that they have no spamt left 
to grapple with the problem itsel f”. 


*T want them to think. [ do not insist on logical or clear or critical 
thinking; just plain thinking’’. 
“The student should also learn that mathematics 1s not completely devoid 
of common sense (1/327 16 = 2)”. 


“Of course we al] have many superior stu’ents, wiose preparation | eaves 
little to be Jesired”. 


(From here on T am on my own). 

Unquestionably, we have a great many students at all levels of 
schooling who exhibit the annoying characteristics and practices of the 
ignorant, the dull, the msfits, the unfits; so that often we feel as 
Mark Twain did about certain people to whom he referred as: ‘They 
are not impossible; they ere highly improbable”. Our common problem as 
I understand it, is t> reduce both in numbers and intensity these so 
disturbing characten stics aid practices. One obvious way - too obvious - 
is to plow under the possessors of these characteristics and practices. 
There are quite a few teachers and administrators wlio regard “ liquida- 
tion” as the only solution. That’s a dangerous doctrine, since liquida- 
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tion knows no bounds. It may begin wth pupils but it may spread to 
teachers and administrators, No, these pupils are here to stay. 

Also, I am convinced, it was no solution for teachers at one level 
to snipe at teachers at anotl.er level. There is enough poor teachiny 
to go arourd, and who of us is free from some poor teaching! I may be 
overly optamistic, but from much of wiat T hear and read T sense a 
growing realization on the part of teachers at all levels that our 
job as very muck a joint enterprise. Let me paraphrase a variously and 
more and more frequently expressed injunction to college teachers by 
college teachers: We must accept our students as we find them. Petter we 
teach them whatever it 1s they failed to learn than spend an equal 
anount of tame and physical and mental energy on complaining alcut their 
high school and grade school teachers and teaching. — 

In conclusion, T want to attempt an explanation anc offer a suygges- 
tion. I know that if T were to ask all the teachers of college mathe- 
matics for just one bit of advice to teachers of high school mathematics 
it would be: “‘Teach for understanding!” JT know that if T were to ask 
all the teachers of higi: school mathematics wiat single factor do they 
stress most in thear teaching, it would be: “un ierstan:ing’’. How shall 
we resolve this paradox? 

In the kand of learnings which requive understanding, there is no 
way to bypass understanding. But - usiny terminology familiar to all 
of us — while understanding thus becomes a necessery condition to 
learning, it is not a sufficient con ‘ition. To put it differently, 
unlerstanding left unattended by drill or practice or application withers 
away. In this sense understanding 1s no different from memorization. 
Sonehow our psychologists failed to impress us with this fact. Perhaps 
they do not know it. Perhaps they “taught” us and we have forgotten. 

Be that as it may - and this is the “punch line” - tt ts quite 
possible that so many students give so little evidence of understanding 
not because they never understood but rather because understandings, 
unlike old soldiers, do die. 

T realize that even the complete acceptance of this exp] anatior 
does not solve our problem, but J believe it points the way. 

Mur best students, and I find substantial agreement on this, —- are 
as good as ever. They understand, they rememler, they keep on leamuing, 
iley are our pride ard joy. (ur poorest students are probably as poor 
as ever. I suppose we’l] fin! a better way of dealing wth them about 
the time pl.ysicians diseever a cure for the “common cold”. Put it is 
that large and ever increasing number of students in “the middle” who 
are really our protlem chal) dre. 

I suggest, I recowaend, that for this numerically dominant part of 
our pre-colle« and college population we reduce the number of topics 
in any given subject. (Here the reference is to’mathematics, but it 
could well mean “any given subject’); that the remaining topics be 
given a fuller, more leisurely, more meaningful] treatment; that a 
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continuity of topics be established so that no topic is ever “ finish- 
ed”; that every topic be made as rigorous as the training and the 
background of the pupil permit: no more an’ no less; that, therefore, 
the transition from topic to topic, from course to course, and even 
from high school mathematics to college mathematics (watever that 
may mean) may be shock proof because of the smoothness of the road 
traversed rather than the use of fancy shock alsorbers. 

It is wth some trepidation that I am about to highlight the above 
suggestions by a reference to a muck discussed document, - the proposed 
six year syllabus.* It 1s a good syllabus; perhaps the best ever proposed 
in the State of New York. I approve of it, as a source book for our 
best students. As a syllabus of minimal essentials it insists on a great 
deal more material than can be taught meaningfully to “ that large and 
ever increasing number of students in the middle’. So long as we have 
State Examinations - Regents, or their equivalent —- the state syl]lebus 
sets the pattern an! the pace. The pattern of the new syl!atus points 
in the right direction; the pace, I fear, wll not relieve the teachers 
from the task of coaching for examinations. And the examinations wil] 
continue to keep the wolf from the doors of the publis! «rs of Regents 
Feview books. 

T won Zer how much and what kind of new evidence do we need to convince 
us that a little leamed well is a great deal more useful to the leamer 
than a lot leamed (7?) poorly; that “covering ground” is often gratuitous 
exercise for the teacher and nothing more. ** 

Give our secondary schoo] teacners half a chance to meet their 
colleagues at the college level half way. “ore than ever the teachers 
of colleyve mathematics are willing to accept the responsibility for 
building the superstructure, provided teachers of elementury and 
secondary school mathematics lay the proper foundation. “Mere numbers 
of stones or brichs. loosely held together do not fom a proper founda- 
tion. 


*“The Six-Year High Schoo} Mithematics Syllabus for New York State.’ 


**Amen Fd. 


Al fred University 
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MISCELLANEOUS NOTES 
Edited by 


Charles K. Robbins 


Artic! -. intended for this department should be sent to Gharles K. Robbins, 
Department of Mathematics, Purdue University, Laiayette, Indiana. 


MAXIMUM AREA IN A CORNER 


"he following theorem seems similar te several isoj<rimetric theorems 
but doesn’t seeu to be readily found in the literature. It furnishes 
the answers to such questions as: 

What 1s the configuration of a simple polygon of maxintm area subject 
to the conditions that two sides meet at a fixed anzle and the other 
n —~2 sides are prescribed in length only? 

T submit it that this apparent gap be fillec. The approach is ele- 

mentary. 
THEOREM: Among all the curves of ;iven length whose endpoints A and B 
lie on two half-lines OL, and OL. inclined at an angle a (by convertion 
the areas included in the hal bewdene swept out by OL, in rotating 
through a clockwise angle c are to be considere? positive.), the one 
which, together with the half-lines OL, and OL,, bounds the maximum 
area (assuming such a maximumizing curve exists) ts a circular are 
whose center of curvature is atOifOQO<a<n7, If 7 Sa < 2n, the 
curve ts a semicircle, 


Lemma [: The AAOB having the maximum area subject to the conditions 
1E = m and angle AOB = ais the one in which OA = OR. 

For of all tmangles having a common base an: equal oppos:te angles, 
the isosceles triangle has the maxinum altitude. 


Lemma T[I: Anony all the curves of given length s whose endpoints A and B 
lie at the extremties of a fixed line segment, the one which, together 
with chord AB, bounds the maximum area is a circular arc. 

See e.g., Polza, “Lectures in the Calculus of Variations, ex. 13, 
p 210, 
Proof of the Theorem: If O<a¢n, the following argument holds. 

Suppose there exists a maximizing curve C of length s. Then C mst 
be a circular arc, since otherwise the area cou]¢@ be increased by 
Lemma JT, and the AAOB wist be isosceles, since otherwise the area 
of AAOP could be increased by Lenma T. Thus C is a circular are with 
its center on the bisector of angle AOR. We now determine the radius r 
so that the area will be a maxinum. 

The area enclosed by the hal f-lines OL, and OL, (meeting at an 
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angle a) and a circle whose center lies on the lasector of a is given by 


P = Sr/2 - r*[sin(S/2r) sin(a/2 - S/2r)]/[sin(a/2) ] 


where S is the length of arc C, and ris the radius of curvature of C. 
For fixed S anda, P is solely a function of r. A necessary condition 
for a relative maximum is that the derivative of P with respect to r 
must vanish. 

Since 


dp Ss 
2 


[sin(a/2 - S/2r)/(sina/2)][S cos S/2r - 2r sin S/2r), 


the derivative vanishes when 
(v) a/2 —~ S/2r = tKn or (ti) tan s/2r = S/2r. 
The area in (1) as given by 
P= S*/(2a + 4K7) for K= 9, t+1+2+43... 


Compatible with positive areas, the abso]ute maxinum occurs when K = 0). 
Consider the function 
S2 
B = 1 [cos( a /2) + cos(a/2 - S/r)]/{sin(a/2)] r>0 


wiich coincides with the values of P whenever (it) is satisfied. It is 
not hard to see that 


< Ss? 9 s? < s° 
B= g—[cos(a/2) + 1}/sin(a/2) = q cot(a/4) — 


by comparing with the series expansion for tan a/4 = a/4+ ... > a/4 
for 27 >a>QO0. 


Since the absolute maxinum of Bis less than the absolute maximum 
of the area given under (i) and the values of P when (it) is satisfied 
coincide with those of B. Then al! the relative maximums of P under 
(it) are less than the absolute maximum under (1). The endpoint maximum 
occurs when r is infinite which occurs as a special case of (it) is 
included therein. The endpoint minimum occurs when r = 0. Therefore 
the absolute maximum occurs when 

S2 
Pp = -—- 


2a 


whence the center of curvature of Cis at O the intersection of the 


two half-lines OL, and OL,. If a $ m™, the above reasoning is valid. 
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lfm <a< 27, then once the area between the solution C and the 
chord ABis fixed the maximum occurs wren AAOB has area zero. This 
reduces the problem to the famous one of the bul}bide and the area 
enclosed between it and the rmver. *ts well-known solution is a semi- 
circle. 

Summary. I[f 0 <a<w, the solution is unique. [f @ = 7, there 
is an infinite nusiter of solutions. Tf 7 < a < 27, there are two 
solutions. 


Vern Hoggatt 


ELLIPSE CONSTRUCTION SHORTCUTS 


P RCBLEM: 
To find shortauts in the process of onstructirg ellipses on pictorial 
lrawings from known equal conjugate ch ameters. 


x 














<_ UF _| — 
FIG. 1 


The basic construction is as show in Fig. 1. Since the diameters 
are equal the major and minor axes will bisect the angle between the 
conjugate diameters and the only problem remaining is to determine 
thear magnitudes. The common construction’ shown accomplishes this in 


that the major axis equals al + bl while the minor equals al - 6l, 
Ref. Engineering Drawing, T. E. French, YcGraw-Hil!; similar con- 
struction may be found in neal Constructive Geometry, Dowsett, Oxford 


Univ. Press and Technical Drawing, Giesecke, Mitchell and Spencer, 
Macmillan. 
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The parr of equal conjugate diaveters 12 and 34, axes xx’ and yy’, 


anc the lanes ab and al, in Fig. 3. 


TO PROVE: 
al + bl al - bl 


eg * -———— and ee * = -- 


CONSTRUCTION: 


Extend line al to 6'(lb’ = 1b), draw b'bf and akf. Construct fd and 


om perpendicular to al. 


PROOF: 


First we prove that the line al will make equal angles (45°) with 


the major and minor axes, see Fig. 2. 


Angles y = 0, since the sides of the angles are perpendicular each 


to each and both angles are acute. 
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+ i= Y and v + 4) = I 


Since the triangle ola is isosceles the angles U and & are equal and 


D+ w+ 90° = 180°, Db + = 90°, 
90)° 

t s Wwe -— gs 45° 
2 


Angle alb = 90° since it 1s insermbed in a semicircle. Therefore both bl 
and b'l make 45° wth bb’ hence bb’ is parallel to yy’. 

Angle bfa is insembed in a semicircle and is therefore 90°. b'f = af 
(opposite sides of an isosceles triangle), therefore the line fd, 
constructed perpendicular to al, bisects ab’. dis not assumed to be 
on yy’. 

fdk and adk are congruent right triangles since af is perpendicul ar 
to yy’, also angles at a and f are 45°, therefore dk is perpendicular 
to af. Therefore dk mst lie along yy’. 

Construct om perpendicular to al, then om = dm.= em = 61/2, and am = 








al/2, therefore am + dm = ad = al/2 + bl/2 = 2  * 9 and am - em = 
al - bl . 
ae= al 2 sm bl /9 Ss Ye . 
2 
APPLICATION: 
Since al + bl and al - bl are, by common geometrical construction, 
) the magnitudes of the major and minor axes respectively of the ellipse 
j |- 2-3-4, then ad and ae represent tle magnitudes of tle semi-major and 


seml-minor axes. 

A line may be drawn thre point 2 parallel to al (45° wth yy’) and 
the distances 2r and 2s measured alon» thas line will be exactly tile 
same as ae and ad, Therefore all other construction may be dispensed 
with other than this 45° line. 

These measurements may then be laid off along axes yy’ and xx’, and 
the ellipses constructed by any convenient method, 


S. B. Elrod 








CURRENT PAPERS AND BOOKS 


Edited by 
H. V. Craig 


This department will present comments on papers previously published in the 
MATHEMATICS MAGAZINE, lists of new books, and book reviews. 

In order that errors may be corrected, results extended, and interesting 
aspects further illuminated, comments on published papers in all departmen's 
are invited. 


Communications intended for this department should be sent in duplicate 


to H. V. Craig, Department of Applied Mathematics, 'niversity of Texas, Austin 
12, Texas. 


Intermediate Algebra, By Paul K. Rees and Fred W. Spar’s. McGraw-Hil |! 
Rook Company, New York, 195]. viaia + 322 pages. $3.25. 


The authors of this book tel] us ir the preface that it is intended 
for a tlree-semester-hour course for students who have had only one 
year of high school algebra. [t 1s, therefore, informal in style and 
offers a great variety of illustrative material that shoul cl ari fy 
the discussions of abstract or technical problems. 

Several features of the book should appeal to teachers and students 
alike: 

1). Problems hav- been gra‘ed and arranged in groups of four of about 
the same difficulty. Thais rnekes it possible to assign problems of 
varying complexity by taking every fourth one, or to concentrate on 
either the simpler ones or the more involved ones as occasion demands. 
Toward the end of each exercise are more difficult problems that ml] 
test the ability of the better students. 

2). There are scattered throughout the book numerous notes that 
point out bits of interesting inforaation (pp. 66, 218, 222, 233), things 
the student should not do (p. 13), little tricks to help him to avoid 
errors (pp. 253, 254, 271), how to recognize the easier of two methods 
of solution (p. 211). 

3). Answers are yiven for the even-numlered problens and for some 
of the odd-numbered problems so that it 1s possible for the student 
to check his om comprehension and accuracy, an! thereby gain con fidence: 
in his ability. 

4). The solution of stated problems is, to many students, the rost 
difficult part of a course in intermediete algebra. The authors have 
Cone an unusually fine jo! in handling problems of this type. They 
not only explain the general procedure tut list specific steps to be 
followed in analyzing the problem and setting up the necessary equation. 
In addition, the examples are divived into steps clearly numbered 
to follow the suggestions already listed (pp. 18-22). The examples 
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are of many different types inclucing “work” , “motion” and “mixture” 
problems (pp. 88-91). The authors make use of the same kind of step -by- 
step explanation in stetement problens leading to systems of linear 
equations (pp. 131-135), ain problems leading to quadratic equations 
(pp. 195-197), and in “ variation” problems (pp. 243-246). Their ex- 
planations and illustrative examples should help the student to master 
the process of translating stated relationslips into correspondiny 
mathematica! symbols. 

5). Purely mechanical problems are presented mt the same detailed 
exp) anations, witl frequent reference to the proper law or rule that 
applies in each step of the solution. The explicit directions given 
by the authors make this an easily understocd textbook. 

The first nine chapters are }uilt aroun’ the equation as the central 
idea. The discussic:: leads from tie formulas of arithmetic, already 
familiar to the student (chapter [), through the fundamental operations 
of algebra (chapter TI), special products and factoring (chapter ITT), 
to fractions and fractional equations in chapter TV. Next the meaning 
of function, functional notation an: graphical representation of 
functions are presented (chapter V) as preparation for simu! taneous 
linear equations (chapter VI). Exponents ani radicals (chapter VII) 
are discussed before the introduction of quadratic equations (chapter 
VIII) and the solution of systems of quadratic equations in two variables 
(chapter IX). In this chapter (IX) special methods have been avoided 
and the method of substitu tion has been emphasized. 

Chapter X deals with ratio, prcportion, and variation; chapter XT, 
with logarithms; chapter XII, with progressions; and chapter XIII, 
with the binomial theorem. The appendix of the textbook includes a 
four-pl ace table of common logarithms, ad a table of powers and roots. 
The tables add to the usefulness of this well written book. 


Vivian Gumno 


Advanced Engineering Mathematics. By C. R. Wylie, Jr. McGraw-Hill, 
N. Y., 1952. $7.50. 


Of beoks on “advanced calcilus” there 1s no end. One wonders, when 
seeing another one, whether its existence is really justified. The 
reviewer is inclined to feel that new books on advanced calcovlus are 
justified because the increasing importance of the subject in the 
cngineering scene makes it necessary to reach a wider circle of engi- 
neering students. Forty years ago when E. B. Wilson wrote one of the 
first adequate American advanced calculus texts the emphasis was on the 
needs of theoretical physicists. Tt was just these physicists who 
became the engineering analysts and demonstrated the importance of 
advanced mathematics to engineering analysis. Today real theoretical 
studies are every day practice in electronics, aeronautics, and other 


branches of engineering. A working knowledge of advanced calculus 
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1s necessary for the comprehension of many important technical papers. 
Pooks on advanced engineering mathematics are necessary. There should 
be a variety of them to meet a diversity of needs. It appears to be thie 
task of a reviewer to assess these needs and the adequacy of the book 
for meeting them. 

I think that ordinary and partial differential equations, Fourier 
analysis, Laplace transforms, functions of a complex variuble, vector 
and tensor analysis, including analytic and differential geometry, 
probability theory, and modern algebraic methods including matrices 
are all necessary for an understandiny of present day engineering 
analysis. [ think that the important part of these subjects can be 
presented within the compass of an 800 page book. The present book 
has nothing on algebra and matrices, probahility theory, tensor malysis, 
or differential geometry. The book does treat methods of numerical! 
analysis, a topic which is better taken up as a separate course along 
with modern digital and analogue computing equipment patterned after 
the delightful little monograph of Hartree and emphasizing the under- 
lying principles. An “ advanced mathematics” book should emphasize 
principles and illustrate with worked out examples. 

What the book does cover it covers well. The emp! asis on problem 
solving is excellent and the worked out examples are instructive. 
As far as teachability is concerned the reviewer cannot comment. His 
viewpoint is entirely that of an applied mathematician working in the 
engineering field. From that standpoint he does feel that the books 
by Burington and Torrance, Jeffreys and Jeffreys, Margenau and Murphy, 
Schelkunoff, and Gilleman should also be used for reference purposes 
although none of these cover probability theory and none of them include 
all the topics which should be included. A book on advanced engineering 
mathematics should include many references for further reading. This 
omission is a common fault and+a bad one because it is impossible and 
undesirable to include al] the detail in one volume. An advanced en- 
gineering mathematics text should be a starting point not a terminal 
point in an engineers study of a mathematica] technique he intends 
to apply in his work. He should be able to use such a book as a gateway 
to a thorough knowledge of the part of the subject he intends to use. 


N. G. Parke 
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PROBLEMS AND QUESTIONS 


Edited by 
C. W. Trigg, Los Angeles City Col lege 


Readers of this department are invited to submit for solution problems 
believed to be new and subject-matter questions that may arise in study, 
in research, or in extra-academic situations. Proposals should be accompanied 
by solutions, when available, and by such information as will assist the editor. 
Ordinarily, problems in well-known textbooks should not be submitted. 

Solutions should be submitted on separate, signed sheets. Figures shoul 
be drawn in India ink and twice the size desired for reproduction. Readers 
are invited to offer heuristic discussions in addition to formal solutions. 

Send «all communications for this department to C. W. Trigg, Los Angeles 
City College, 855 N. Vermont Ave., Los Angeles 29, California. 


PROPOSALS 


147. Proposed by Leon Bankoff, Los Angeles, California, 


In a triangle with integer sides, the side opposite the 120 angle 
is 1729. Find all possible values of the pair of other sides. 


148. Proposed by D. L. MacKay, Manchester Depot, Vermont. 


Upon the sides of triangle ABC the squares ABDE, BCFG, ACHL are 
constructed exterior to the triangle. Construct triangle ABC given the 
points A’, B', C' which are the intersections of DE ad HL, ED ad FG, 
GF and LH, respectively. 


149. Proposed by L. C. Lay, John Muir College, Pasadena, Calif. 


A has a gambling device so arranged that he always wins 3 times 
out of every sequence of 5 plays. B suspects that the game is “ fixed” 
and demands that A always wager half of his resources, which B will 
match with an equal amount. B’s funds are assumed to be unlimited. 

(1) Show that A always loses over a sequence, regardless of the 
order of his 3 winning and 2 losing plays. 

(2) Determine the best fixed percentage of his resources for A to 
wager if he wins 3 out of every 5 plays. 

(3) If A mst always wager 1/2 of his resources, determine to the 
nearest integer the number of times he must win out of a sequence of 
one hundred plays in order to break even. 


150. Proposed by P. D. Thomas, U. S. Coast and Geodetic Survey. 


The lines joining the vertices of a triangle to the intemal points 
of contact of the escribed circles meet in a point. The perpendicul ars 
upon the sides of the triangle from the excenters meet in a point. Show 
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that these two points together with the orthocenter and the incenter of 
the triangle are the vertices of a parallelogram. 


151. Proposed by Dewey Duncan, East Los Angeles Junior College. 

In a recent text on Backgrounds for Secondary Mathematics Teachers 
the following statement appears: ‘‘Consider the sequence 1, 2, 3, 2 %, 
2, 1%, 1%, 2, 2%, 2 1/8, ***. It has the limst 2, as can be shown 
by locating the values on a number scale. There is no single formula 
for this sequence.” Refute this last assertion by example. 


152. Proposed by Malcolm Robertson, Rutgers University. 


In fin@ing the area, A = 7ab, of the ellipse 0? = a*b*/(a*sin*? 
+ b*cos*4) a student gets an incorrect answer by proceeding as follows: 





en ab? ("a sec’9 d9 
alts x | niakis 2 | b2 + (a tan 9)? 
0 0 
ab2r] 2n abr 
= ar arc tan (; tan | = —l!larc tan 0 - arc tan o = Q. 
2 [b . 2 | 


Detect and explain the source of error. 


153. Proposed by John R. Hatcher, Brown University. 


Prove (without using the exponential function ef'?)) that if f(z) = 
u+ iv is entire and u# 0, then f(z) is a wnstant. 


SOLUTIONS 
Sums of Sets of Binomial Coefficients 


117. [Nov. 1951] Proposed by H. D. Grossman, New York, N.Y. 


If n is a positive integer, prove that: a) The maximum di fference 


n n n n n 
between any two of the four sums, | + \+ t+ oee) ("| " 
ay, 


0 4 8 


x « Cnr » Sb bo eae Oe se 


where [n/2] means the greatest integer < n/2. b) The maximum di ffer- 


n n n n 
ence between any two of the six sums, | + | | + | | + +. K 
0 4 12 8 


+ ("| 2 | ae, eae f F sal ; (| i ie Bares gin/2)-1 if 
7 13 5 li 17 


gforti 7s ig n is odd, 


nis even, or 
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Solution by E, P. Starke, Rutgers University. a) Let S,, S;, 
S,, Sz, respectively represent the four given sums. From the binomial 
expansion of (1 + x)" with x = it, -l, -i, 1, we have 


dtd td) - Gy 
, 4 
J 7h Tt 
ol Gets Ge 


Upon multiplying these equations, respectively, by 1%, i*, i>*, 1 
(a = 1,2, 3, 4) and adding, we obtain 


(1 + i)" = 


to 


(1 - i)" = 


| 
aeare| 
| 


! 
’ a 
— | 
UH eS 
+ 
| 
ae 
! 


w > 


(3 + oF 


+ 


" 


i(l+i"- (1- i) ® + = 4S., -(1+ i)" -(1- i)®+ 2" = 


! 
— 
” 


3’ 


“i(L+ip® + a(l-i)* + 24S, (Leip +(1-i)* + P= 4S 


0’ 





S;, - S.|, we find, after obvious sinpli fi- 


A) l 4 D = 
whereupon, putliing ij j 


cations, that 


&, «2,4 ogy of if life. @ 7 {Ql +i)" +(1-1)"|/2, 


O11 0 
g= [CL + iv??? + (1- i)"*?]/4, Di g=D,,= [CL + iv) + (L- i?" '|/2. 


The possible values of the D’s are now easily computed. For n = 2r, 
they are 0, 2""', 2"; for n= 2r +1, they are 0, 2 only. From this, 
the desired result follows directly. 

b) In a analogous fashion we set down the expansions for (1 + w& 
(ec = 1, 2, ***, 46), where # = cos(7/3) + i sin(7/3). If these are 
multiplied by w%, uf -++ > aP® = 1] and added, there results 


ry. 


w*(] + w)"+ we *(] +w*)" + 4%] + w*)" + y%(] +w)" . 9” ie 65. 


for a= 1, 2, ***, 6. Expressions for the fi fteen D; ’s can now be 
formed and will be found to be of three types: 


1) 6D), = |2(1 + @)" + 21 + w)"!, 


6D, = | 2e2(1 + w)™ + 204(1 + w5)"|, 





= | 21 +a)" + Qa4(1 + 5)"]; 
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(2) 6D, » =|aP(1 + w)" + w3(] + po)h o Mtl 4 pft2 4 5nt4 $ yot5| 
6D, 5 = |? (1 + w)™ + aA (1 tw) e PTS 4 whe PTl 4 5M) 


D D 


34’ 45’ D 


and analogous results for D bai 


05’ 


(3) 6D,,= [v2 (1 + wy PPT + aS (1 + 5)Bth te oe atl + 544 SMFS 


35’ Dy 4» D 


- #@ +] = 0 are useful in simplifying.) 


and analogous results for D,,, D 
2 


15» Dyo- (The relations 
w = -j, w 

It will save space to put H, for w* + w°*, Tf n= 2r, the D's take 
the following forms: 


- ~f = - ‘ =e : 
a Oo oes i st ee Og hs 
6D, 2° 13H, 45 * Hares *Harsgls SDa5* 37H, 42 * Horas tHe, le °° * 
6D,, > |3°(H,,. * Hpy3) * He, * Harayl, 


For n = 2r + 1 we have: 


D,,=3'"'|H,+H,,,|, D,,=3”'|H,,,¢H D,.= 3" '|H,.,,+H 


yeyle 25 PF: 


6D,, = |3'(H,,, * H.44) + H +H 


r+4 2rt2 aregls °° 


6D,. = |3°°'H,,, + H +H 


2r+) areal 


Since there are only six distinct powers of # it is easy to see 
that, whatever the value of a, H, must be +1] or +2. Also H, * Has 
can be only +3 or 0. With these facts in mind and remembering that 
the D’s are all integers, we find the only possible values of the D's. 
They are, 


for n= 2r: 2°3""', 37°', (37 + 1)/2, (37°' + 1)/2, 0; 
for n= 2r +1: 37, (3" + 1)/2, 0. 


In every case the maximum value is actually attained by one of D,,, 


D D 


14? 25° 


If these results are assumed at the outset, it is not difficult 
to establish them by induction on n. However, the form of Ds; de- 
pends on the residue of n mod 12, so that induction requires working 
from n to n+ 1, n+ 2, ***, n+ 12 before it is complete. 

If the binomial coefficients of order n are grouped into five sums 
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1) a manner analogous to tlie above, the maximum D; equals the (n + 1)st 
Fibonacci number (1, 1, 2, 3, 5, 8, °°; Snes = f, + Sno) [See 
National Mathematics Magazine, 13, 292, (Varch 1939)]. 

Other references dealing with the sums of sets of the binomial 
coefficients are: American Vathematical Vonthly, 39, 304, (1932); 
45, 320, (1938); and National “Mathematics Vagazine, 10, 145, (February 
1934). 


Equation Leading to a Fermat Equation 


121. (January 1952] Proposed by Norman Anning, University of Michigan, 


Solve in positive integers, (x + iy)? = x + (a pure imaginary). 
For instance, (7 + 4i)° = 7 + 524i. 


Solution by H. R. Leifer, Pittsburgh, Pa. Let (x + iy)? = x + im, 
then x° - xy? + i(3 x*y-y?) = x + im, Equating the real parts, we have 
x? - 3xy" = x, whereupon eit er: 

(1) x = 0, which leads to the trivial non-negative solutions, 
y equals any positive integer; or 

(2) x* - 3y? = 1. This is a Pell (more properly, Fermat) equation, 
which has an infinity of solutions. Consider 


ae ae 


JS a | O ame cee oe oe * 


where the successive convergents are 


pz 3 TW SB THT MS |US 

1’ 1’ 3’ 4° 11" 15" 41° 56° 153’ 209° 
The numerator of each even convergent gives a value for x and its 
denominator the correspo:dirg value for y. Thus the first five solutions 
are (2,1), (7,4), (25,15), (97,5), (362,209). Or, the general solution 
1s given by 


x= [(2 +73)" +(2-/3)")/2, y= [(2+73)" - (2- 73)")/273, 
ae (@ £,; & *°*}. 


Also solved by Leon Bankoff, Los Angeles, Calif; A. L. Epstein, 
Cambridge, Mass.; M. S. Klamkin, Polytechnic Institute of Brooklyn; 
F. L. Miksa, Aurora, Ill.; and the proposer. 

Anning gave the recurrence relations x, = 2x%,_,; + 3y¥,-1, Y¥, * 
X,-; * 2y,-, starting wit! the obvious solution (1,0). He also noted 
that u, = 4u,_, - u._, applies to so0th the x’s and the y's. Banko f f 


gave x, * X92 ° bY ny Yn * 2xn-4 * Yn-2° 
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A Circular Cubic 


122. [January 1952] Proposed by P. DL. Thomas, U. S. Coast und Geodetic 
Survey, Washington, D.C. 


Show that the envelope of a circle, the square of whose tangent 
from tle origin is equal to the ratio of the abscissa to the ordinate 
of its center, the center lying on the parabola x = ay*, is a cireilar 
cubic, one of whose asymptotes is parallel to the y-axis. 


Soluticn by R. E. Jackson, Student, University of California at Los 
Angeles. The circle with center at (6,c) has a radius equal to 
(b2 + c2 - t?)*% where the length of the tangent, t, from the origin 
is given by t* = b/c. Hence the equation of the circle is 


(x - b)* + (y - c)? = b? + c® - b/e. 
Since b = ac’, the equation of the circle assumes the form 


f(c) = x2 - Qac*x + y? - Ay + ac = O. 
Then 
df/dce = -4acx - 2y+a=Q0. 


The equation of the envelope of the family of circles is ol)tained by 
eliminating c between f(c) and df/dc. Thus 


Rax(x? + y*) + (2y - a)? = O. 


The terms of the third degree equated to zero give the lines through the 
Origin parallel to the asymptotes, that is, x = 9 and x = tiy. Hence 
the cubic is circular. Writing the equation of the cbic in the form 
y*(Bax + 4) - day + Bax? + a* = 0, we see that the asymptote par} lel to 


the y-axis is x = -]/2a, 


Also solved by M. S. Klamkin, Polytechnic Institute of Brooklyn; 
A. Sisk, Maryville, Tern.; and the proposer. 


Altitude of a Spherical Triangle 


123. [January 1952] Proposed by Joseph Barnett, Jr., Clarksturg, W.Va. 


Theorem: A necessary and sufficient condition that the foot of a 
perpendicular from the vertex of a spherical triangle to the circle 
containing the opposite side fal] or that side is that the angles 
adjacent to that side be of the same species. 


Solution by the Proposer. Let the side BB’ of the spherical tri- 
angle BB'A be contained by the great circle DBD'B’, with pole P, 
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and with D (between B’ and B) and D' (between 6 and B') being the inter- 
sections of this great circle with the great circle through P and A. Then 
AD and AD’ are the only possille perpendiculars from A to the containing 
circle. 

In triangle BAP, BA+ PA > PB = PD = PA+ AD, so BA > AD. Then 
90° = angle BDA > angle ABD. Simlarly, 9° > angle AB'D. It follows 
that angles ABD’ and AB'D are obtuse. This proves the necessary con- 
dition. The converse is true, for if the adjacent angles are of the 
same species, the perpendicular must lie wthin the triangle. 


Composite Resistances 


124. [Jauary 1952] Proposed by Leo Moser, University of Alberta, Canada. 

Prove that if p and q are positive integers not exceeding the integer 
n, then 1t 1s possible to arrange n or fewer unit resistances to give 
a combined resistance of p/q. 


Solution by the Proposer. We proceed by induction. The theorem is 
clearly true for n = 1 and also for p = q. 

If p 7 q, one of the fractions p/(q -— p) or (p - q)/q will have 
both numerator and denominator positive and less than n. “ence one of 
these will be obtainable using n - 1] or fewer unit resistances. (This 
is the induction hypothesis.) If the first of these 1s obtainable, then 
we get a resistance of p/q by putting it in parallel wth a single unit 
resistance, while if the second is obtainable, we get a resistance 
of p/q by placing it in series with a unit resistance. 


Also solve:' by George Baker, Student, California [Institute of 
Techno logy. 


Problem Related to a kicatti kLquation 


125. [January 1952] Proposed by William Leong, University of California 
at Berkeley 


2 

Consider the sequence of numbers {a } where 3a, = 1, 7a, = a, , 
lla, = 2a,a.,, 1Sa, = a, + 2a,a.,, 19a, = 2(a,a, + a,a.), 23a, = 
a, + 2(a,a, + a,a,), °°" (1) 


4 ; ; 
Let u «= lim a /a . Then show that (a) the number U exists, and 
n~@ n+1 


(6) u satisfies the equation 








k 4k 
@ -1) 
5 ~ sa te a = 9 (2) 
k=} 3°4°7°8*** (4k - 5)(4k - 4)(4k - 1) 4k 
Solution by the proposer. Consider the Ricatti equation 
y' = y? + x? (3) 


If we set y = - u'/u then u” + ux* = 0. 
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The latter equation has, by proper choice of constants of integration, 
the series solution 


4 
1 - ¢,x + Cy 


u x - Cc 


a @ 0-8 (4) 


where c, = —, C, = 
3°4 3°4°7°8 








gs : 
mM 3+A*7T Re ++ (4R-5)( 4e- 4) (A4R-1) 4k 


Hence the original equation (3) has a solution 





3 c; - 2c,x4 + 3c ,x° Soe z 
y= 4x re (*) 
4 8 12 
1 -¢,* + cox - C 4x + 
~ 3 7 11 15 
a,x + a,x" + a,x + a,x + R (4) 


where the a; are given by (1). 
[The direct division involved in (5) leads to a different set of 
recurrence formulae fron (1): 


n _ 
(-])"q = -4ne +c & «ft 


h noi" n= 2°72 


e van, of (-1)"c a + (<1) ¢ a ° 
2n-2 lone l 


The formulae (1) are more ccnvenient for calculating the a,; they were 


derived by working with (3) instead of the transformed equation. | 

We now raise the question of the radius of convergence of series 
(6). Ry the ratio test we have that (6) converges for all values of x 
for waich |x| <u, On the other hand, the series (4) and its di ffer- 
entiated series represent entire functions by the ratio test. ‘lence 
it follows from (5) that if u(x) has a zero then U exists and must 
be the smallest positive zero of u(x). Fut it is easy to compute 
u(3) = -0.33 and u(() = +1. Hence the statement follows. 

1 ] 2 13 46 


The first eleven a; are: =’ =’ 


332-7) 33-791) 34e8e7 291] 355729119 





15178 404 190571 
36-5-73-1] 7-19-23 37 -5+77+1] 7-10-93 37 «52-74-11 7-19-93 4] 











5 858 822 ____1 87 447 754 7 
39 +5 3-74-11 3+19+93°3] 310-534-7541] 53+13+192+93°3] 








4389 246 452, — tC 
3'° .59.75.11% 13-107 .23.31-43 


With the aid of a ten-place logarithm table, I have computed the 
first ten ratios a,/a,., and have found that the numbers become stable 





fairly rapidly. In fact, 
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Genoese 4 . - 
Va,/a,, = ¥15 194 854 857/943 723 877 = 2.00314 742 -*-, 
a 4. 








¥a,,/a,, = 740 580 126 711 / 2 520 355 486 = 2.00314 736 °°, 


Also TIT computed u(2.003147) 
u (2.003) 48) 


+0.00000 050°°*° 
-0.90000 090° +:, 


Note that at a zero of u(x) the curve of u(x) vs. * has a point of 
inflection since u” = -ux*, This fact facilitates the location of the 
zero. 

The series (6) was first dernved Ly Jeunes Remoulla in 1703. Reference 
is made to a letter to Leibniz (Leibniz: Gesamellte Werke, 1855, Vol. 3, 


p. 75.) This work does not seem to disauss the radius of convergence. 


Target Practice 


126. [March 1952] Proposed by George Pate, Gordon Military College, 
Georgia, 


In firing a mfle at a target from a given ‘distance, suppose that 
the probability of hitting the bull’ s-eye is 0.3. What is the smallest 
number of shots which must be fired in order that the probalility of 
hitting the bull’ s-eye at least once mll be 0.9? 


Solution by L. D. Rice, Timken Roller Rearing Co., Canton, Ohio. 
The probability of not hitting the target with a single shot is 1 -— 0.3 
or 0.7. Then 0.7" = (1 - 0.9) = 0.1, wiere n is the number of shots 
ani 0.1 as the probability of missing every time. Hence, n log 0.7 
= log 0.1 and n = (-1.9000)/(-0.1549) = 6.46. Tt follows that seven 
shots must be fired, 

Also solved by F. F. Dorsey, South Orange, N. J.; J. M. Howell, Los 
Angeles City College; W. S. Klamkin, Polytechnic Institute of Brooklyn; 
Joel Levy, University of Michigan; C. D. Smith, University of Alabama; 
and J. A. Tierney, U. S. Naval Academy. 


Circle Inscribed in the Arbclos 
127. (March 1952) Proposed by M. R. Watson, San Fernando High School, 
California, 


AB is a dianeter of circle (0). The externally tangent circles 
(O") and (O") have their centers on AR and are also internally tangent 
to (O). Find the radius of (O“) which is tangent to (0), (0'), and (0"). 
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I. Solution by Howard Eve,, Champlain College. Denote the radii 
of circles (0), (0'), (0"), (O") by r, r,, 
reat an’ let & denote the non-obtuse angle formed by AB and OO". 
Then by the law of cosines, 


ro, Ts respectively. Tate 


(r, + r,)° aur « r,)° o tr = r,)° + 2(r —-r,)(r -1r,) cos5, 


2 2 2 
(r, + 3) = (r-r,) + (r,-r.) - 2r-r,)\(r-r,) cos 6. 
Eliminating cos 7, we find 


2 2 2 
l(r, +r,) -(r-r,) -(r- r,) I(r - r,) 


2 2 2 
+((r, +r.) -(r-r,) -(r- r,) ]l(r - r,)= 0. 


3 
Then, using the fact that r, + r, = r, we obtain 
2 2 2 2 2 
. »_ 
r. (r, + r)ryro/(r, +rrier, ) or 2rr,r,/(r Ly o'*, bs 


II. Solution by A. Sisk, Maryville College, Tennes.ce. In the triangle 
C'O"O", using OO” as the cevian, by Stewart’s Theorem we hav: 


(r-r,)*(r, + 7.) = (ry + 75)*(r - ry) + (rg 4 g)*(r = 7) 


- (r, . r,)(r -r,)(r-Tr,). 


+ Pes 


s 2 
Solving for r,, we secure r, = (r, + r,)rir,/(r,° + ryr, 
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I5I. Solution by M.S. Klomkin, Polytechnic Institute of Brooklyn, 
Appl yang Heron’s formula *o triangles 0'0"0", O'00", O"O0" ad equating 
areas, we get 


v(r+ r,)r,(r- 2 di rs) evrr rie-r,- r.) +Vrror.(r- ry" r,). 





Clearing of radicals and using the relationship r= r, + r5, we obtain 


rm 2 r.) - 3r*r.*r.? = O, 


1 2 


2409 2 
r= - rf + Orreor - 
») 2rr, 9r4(r r 


1 


a _ 2 a 2 
Thus rz = rryr,/(r rir.) or ryra(r, + ra)/(r,? + ryry tr 


2 
1 2 ) . 


2 

IV. Solution by A. L. Epstein, Cantridge Kesearch Center, Massa- 
chusetts. Jn a system of coordinates with the origin at the point of 
tangency of (O') and (O") amd the A-axis coinciding with AR, the centers 
of the circles will be O(-r, + r,, 0), O'(-r,, 0), o"(r,, 0), and 


O"(x, y). Then by the distance formula, we have 


19°C")? =:4s.¢ rg) ¢-9* 2 fe, 7)", 
Ow * «42 <9f)* + fF = Gs 457, 
(00")? = (x + i, = Pi + y* = (r, e te 2 Pee 


Now the necessary and sufficient condition tlat these three equations 
in x? + y*, x, and 1 have a non-trivial solution is that the deter- 


tanant of thear coe ficients vanish. That is, 


l 
9 9 ~ 
l 2r (2r, + r,)(-r,) | 
-9 9 “ = 
] 2r, (2r, + r,)(-r,) 0. 
9 ‘a 9 = a 
l Ar, ro) (2r, r,)(-2r, + rz) 
Whe = 2 + 2 " o 
hereupon r, = ryrj(r, + r,)/(r,° + rir, + r,°). There are two tri 


tangent circles with this racius and symmetrical to AB, 


V. Solution by Leon Bankoff, Los Angeles, California. Each of 
the two congruent parts into which AB divides this configuration is 
an arbelos with an inscribed circle. Py a theorem of Pappus [ see 
R. A. Johnson, Modern Geometry, Houghton-\u fflain Go. (1929), page 117} 
the perpendicular O"H from O” upon AB is equal to 2r,. Hence the area 
of triangle 0'0"0" is 2r,(r, + r,)/2. Py Neron’s formula this area 


alsoisVrirjr(r, + r, + rz). When these expressions for the area are 
; 2 
equated, and the equation is solved, we obtain r, = ryrj(r; + rp)/(r,° + 


2 
r + 
"2 * T°). 
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Another approach can be made by using a theorem of Casey [A Sequel 
to Euclid, Dublin (1884), page 118]: If a variable circle toucbes two 
fixed circles, its radius has a constant ratio to the perpendicul ar 
from its center on the radical axis. If we consider (0) and (0”) 
variable, 


r,/AH = r,/AQ" or AH = r,(2r, + r.)/r, . 
If we consider (0’) and (0”) variable, 
r,/HB = r,/BO' or 4HB =r,(2r, + r,)/r, 
Since AH + HB = AB = 2r, we have 
2r = r,(2r, + r,)/r, + r3(2r, + r,)/r, 


| ‘ 2 2 
whereupon r, rrr,/( r,° + ryr, + 7r,°) as before. 


Stall mother general method follows the specific one of Archimedes 
given in OEUVRES D’ARCHIMEDE (literal trans! ation by F. Peyrard), 
Francois Buisson, Paris (1807), p. 431. Let PQ be the diameter of 0 
parallel to AB; R and S and D respectively be the common points of 0",0' 
and 0”,0” and 0',0"; I and L respectively be the intersections of 
CA,(0') and CB,(0"); M and N respectively be the intersections of 
DI,AQ and DL,PB; V and W respectively be the intersections of PM,AB 
and QN,AB. Triangles CO"P and OOA are isosceles and similar so CA 
passes through P. Likewise, CB passes through Y, QA passes through 
R, and PB passes through S. Now angles inscribed in semi-circles are 
right angles so PRD and QSD are straight lines. In triangle APD, 
AR and DI are altitudes, so PV 1s an altitude also. Likewise, QW is 
also perpendicular to AB. It follows from the properties of parallel 
lines that 


AD/DB = AM/MQ- AV/VW or 2r,/2r,= AV/2r,, and 
AD/DB = PN/NB = VW/WB or 2r,/2r, = 2r,/WB. 
Then 


2r 
Finally 


,* @ = AV + VW + WB = 2r,r,/r, + 2r, + 2r,r,/r,- 


r,=(r, + r,)/(r,/r, + 1+ 1,/r,) or (r, + ro)ryro/(r,* * Ff, * ver 


A simple construction of (0") proceeds as follows: Connect C, the 
extremity of the radius O”C perpendicular to AB, and 0’. At D, the 
point of tangency of (0) and (0"), erect a perpendicular to AB meeting 
Co’ in E. On O’O' lay off OF = DE. Draw FE, and O"G parallel to 
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FE with G falling on O'C, Let K be the foot of the perpeidicul ar from 
G to AB. With O' as center and radius equal to O'D + GK, and with 
O° as center and rahius equal to O"D + CK describe circles. With the 
intersection of these circles as center and radius GK describe (("). 

The proof of the construction starts with the similar triangles 
Epo’ and CO"O', in which ED/DO’ = Co"/0"0', so ED= ryri/(r, + r,). 
Then in similar triangles GO"O’ and EFO', GK/O'O” = ED/O'F, so 


GK = ((r, + ro)ryr,/(r, + r,))/lr, + ro- ryrof(r, + r.)] 


or ryro(r, + r,) Kr,* + rir, + E> 
Also solved by Leon Bankoff (in a fourth way using Method ITI), 
and by C. D. Smith, University of Alabama, 


Editorial note: This is a special case of the general problem of 


Apollonius. If three circles with radii r r, touch each other 


- 
eee deed 

externally, then the radii of the two circles which touc'’ the three 
circles are 


icliomasancngeiticninaiit 
_ / + 
r r Pera /lavr ror (r, r 


+ 4 + + 
. Pr.) etre ror r.r,)] 


2 23 3 
contrary to the statement made in E. W. Hobson, A Treatise on Plane 
Trigonometry, Cambridge, Fourth Edition (1918), pa,e 216, problem 2. 
If we use the negative sign and impose the condition that r= r 
we obtain the result of the qrrent problem. 


+r 


1 2” 


QUICKIES 


From time to time this department will putlish problems which may be solved 
by laborious methods, but which with the proper insight may be disposed of 
with dispatch. Readers are urged to submit their favorite problems of this 
type, together with the elegant solution and the source, if known 


@ 15. [Sept. 1950) Find the sum of the squares of the coefficients 
in the expansion of (a + 6)"." M. S. Klamkin offers this alternative 
solution. (1 + x)" = (f) + (P)x + se + (A)#, and (1 + 1/x)” = (9) 
+ ()/x + +++ + (*)/x". Multiplying we have that the sum of the squares 
of the coefficients is the constant term (middle term) of (1 + x)*"/z". 
That is, Qn! /(n!)? or (2M), 


Q 70. (September 1952). The last sentence of A 70, page 54, should 
read: “Therefore the maximum value scught is 


| : a2 |"*"/0N +i a i» 
1 n n=1 n 
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Q@ 71. Prove that the moment of inertia of a square about a line in the 
plane of the square and through its center is independent of the onen- 
tation of the line. [Submitted by Leo Moser. | 


@ 72. Evaluate: lim 2n/csce(7/n). [Submitted by Charles Salkind, | 


n-@® 


Q@ 73. Find the sum of the squares of the roots of 3x4 ~ 36x° + 15x’ 
-~ 1952 = 0. [Submitted by D. E. Thoro.] 


@ 74. In Ripley’ s “New Believe It or Not Book” the following problen 
is given: “Wiat number if divided by 10 leaves a remainder of 9, divided 
by 9 leaves a remainder of 8, divided by 8 leaves a remainder of 7, 
*++) divided by 2 leaves a remainder of 1.” The answer given in that 
book is 14,622,047,999. Find a smaller solution. [Submitted by Lec 
Moser. |] 
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